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Abstract 

We present a brief overview of the most relevant current issues related to quarkonium production in high energy 
proton-proton and proton-nucleus collisions along with some perspectives. After reviewing recent experimental and 
theoretical results on quarkonium production in pp and pA collisions, we discuss the emerging field of polarisation 
studies. Thereafter, we report on issues related to heavy-quark production, both in pp and pA collisions, comple- 
mented by AA collisions. To put the work in a broader perspective, we emphasize the need for new observables to 
investigate quarkonium production mechanisms and reiterate the qualities that make quarkonia a unique tool for many 
investigations in particle and nuclear physics. 
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1. Introduction 

The attention devoted to heavy quarkonium states 
started with the discovery of the Jlif/ charmonium me- 
son in 1974 followed by the T bottomonium meson 
in 1977. From the theoretical point of view, quarko- 
nium bound states offer a solid ground to probe Quan- 
tum Chromodynamics (QCD), due to the high scale pro- 
vided by the large mass of the heavy quarks. 

The application of perturbative techniques together 
with the factorization principle gave birth to the Color- 
Singlet Model. Since then, the appearance of puzzling 
measurements has never stopped and has led to new 
challenges for theorists. These puzzles required the in- 
troduction of new ideas providing new probes for the 
understanding of QCD. Fifteen years ago, the observa- 
tion of an excess in charmonium production reported by 
the CDF Collaboration, by orders of magnitude over the 
theoretical predictions available at that time, gave rise to 
the theory of Non-Relativistic QCD. 

Data collected at the Tevatron, at HERA, and at low 
energy e^e^ colliders has never ceased to challenge 
the existing theoretical models: the apparent violation 
of universality arising when comparing data from the 
hadron-hadron and the lepton-hadron colliders, the dis- 
agreement between the predictions for the polarization 
of the J I if/ produced in hadronic collisions and the cur- 
rent data, as well as the excess of double charmonium 
production first observed by Belle. The solution to these 
puzzles requires new theoretical developments on the 
underlying production mechanism(s), as the computa- 
tion of higher-order corrections for characteristic pro- 
cesses and the study of new production processes, in 
addition to investigations of new observables not ana- 
lyzed so far 



The interest in this field and its developments are not 
limited to the issue of the production mechanisms. The 
progress in lattice calculations and effective field theo- 
ries have converted quarkonium physics into a power- 
ful tool to measure the mass of the heavy quarks and 
the strength of the QCD coupling. The properties of 
production and absorption of quarkonium in a nuclear 
medium provide quantitative inputs for the study of 
QCD at high density and temperature. In fact, charmo- 
nium production off nuclei is one of the most promising 
probes for studying properties of matter created in ul- 
trarelativistic heavy-ion collisions. Since quarkonium is 
heavy, it can be used as a probe of the properties of the 
medium created in these collisions, such as the intensity 
of interactions and possible thermalization. 

Lattice QCD calculations predict that, at sufficiently 
large energy densities, hadronic matter undergoes a 
phase transition to a plasma of deconfined quarks and 
gluons. Since 1986, substantial efforts have been ded- 
icated to the research of high-energy heavy-ion colli- 
sions in order to reveal the existence of this phase tran- 
sition and to analyze the properties of strongly interact- 
ing matter in the new phase. The study of quarkonium 
production and suppression is among the most interest- 
ing investigations in this field since calculations indi- 
cate that the QCD binding potential is screened in the 
QGP phase. The level of screening increases with the 
energy density/temperature of the system. Given the ex- 
istence of several quarkonium states, each of them with 
different binding energies, it is expected that they will 
sequentially melt into open charm or bottom mesons 
above succesive energy density thresholds. Moreover, 
the SPS and RHIC data on charmonium physics have 
brought to Ught further interesting results, among them 
puzzling features in proton-nucleus data, which high- 
lighted new aspects of charmonium physics in nuclear 
reactions, namely the role of cold nuclear matter effects. 

The startup of the LHC and the opening of the new 
energy frontier will offer new and challenging possibil- 
ities for the study of quarkonia. 

All the above reasons have motivated the organization 
of the Quarkonium 2010 workshop held at the Ecole 
Polytechnique (Palaiseau, France) from July 29 to July 
31, 2010. This workshop, gathering both experimen- 
talists and theorists, was devoted to finding answers to 
the numerous quarkonium-hadroproduction puzzles at 
the dawn of the LHC era and the concurrent apogee of 
the Tevatron and RHIC. Introductory and review talks 
focusing upon recent theoretical and experimental re- 
sults were presented, in addition to six topical Working 
Group (WG) discussion sessions, each one devoted to a 
specific issue: 
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• WGl: Quarkonium production in 

• WG2: Quarkonium production in pA 

• WG3: Polarisation and Feeddown 

• WG4: Open Heavy Flavour (vs hidden) 

• WG5: Quarkonium (production) as a tool 

• WG6: New Observables in Quarkonium produc- 
tion 

We present here a comprehensive review addressing 
all these matters. The following sections summarize the 
working group discussions, together with the achieved 
conclusions. In section 8 we prioritize directions for 
ongoing and future developments. 

2. Quarkonium production in pp collisions 

Among the wealth of quarkonium production mea- 
surements, the first CDF analyses of direct J /iff and 
^{2S) production' at = 1.8 TeV [1, 2] are likely 
the most important to date. They revealed that the mea- 
sured rates were more than an order of magnitude larger 
than leading order (LO), color-singlet model (CSM) cal- 
culations [3], believed at that time - the mid 1990's- 
to be the most straightforward application of perturba- 
tive QCD to quarkonium production. Both these dis- 
crepancies motivated a number of theoretical investiga- 
tions on quarkonium hadroproduction, including within 
the NRQCD factorization framework [4]. In the latter 
approach, quarkonium production can also proceed via 
creation of color-octet QQ pairs, present in higher-Fock 
states, whose effects are believed to be suppressed by 
powers of the relative QQ velocity, v. 

Despite advances, there is still no clear picture of the 
quarkonium hadroproduction mechanisms. Such mech- 
anisms would have to explain both the cross section and 
polarization measurements at the Tevatron [1, 2, 5, 6, 7, 
8, 9] and at RHIC [10, 11, 12, 13, 14, 15]. Obviously, 
the mechanisms involved in hadroproduction should 
also comply with the constraints from photo/electro- 
production (see [16] and references therein) and pro- 
duction in e^e^ annihiliation (see [17] and references 
therein). Here we discuss the hadroproduction cross 
section only and leave the discussion of polarization 



'"Prompt production" excludes quarkonium production from 
weak decays of more massive states, such as tlie B meson. "Direct 
production" furtlier excludes quarkonium production from feeddown, 
via the electromagnetic and strong interactions, from more massive 
states, such as higher-mass quarkonium states. 



to Working Group 3 (Sec. 4). We only note that 
the NRQCD approach, based on the dominance of the 
color-octet transition, so far fails to provide a consistent 
description of production and polarization as well as in- 
clusive production in e^e^ annihiliation. This failure 
may be because the charmonium system is too light for 
relativistic effects to be small and thus the velocity ex- 
pansion [4] may have been truncated at too low an order. 
If the convergence of the velocity expansion of NRQCD 
is the issue, then one would expect better agreement be- 
tween theory and data for the T. This could explain why 
only color-singlet contributions [18] (LO in v) appear to 
be in better agreement with Y [19, 20, 21, 22] than ^ 
production [9]. 

At collider energies, quarkonium production occurs 
predominantly through gg channels. Higher order as 
corrections to the S states have only recently been cal- 
culated [23, 24, 25, 26, 27]. The results show that the 
total cross sections do not increase much when these 
corrections [24, 28] are included so that the perturba- 
tive series seems to converge, except perhaps at large 
where initial-state radiation effect may need to be 
resummed [29]. However, the color-singlet corrections 
at high pt are very large because the QCD corrections 
to the CSM open new production channels important at 
high pj. Similar behavior has also been seen in photo- 
production [30]. 

At LO (a,), the cross section differential in scales 
as pj^^ while several different diagrams which contribute 
at NLO decrease as pj^ or pj^. At NLO in color-octet 
production, the higher-order corrections to S state pro- 
duction do not substantially harden the pj distributions. 
The LO color octet channel, gluon fragmentation, al- 
ready scales with the smallest possible power of pj, 
Pj^. However, there are substantial fragmentation con- 
tributions to '5o production which enhance high pj pro- 
duction of these states. At NNLO (a,), important new 
channels also appear. Color singlet gluon fragmenta- 
tion is relevant in the limit pj s> m. Conversely, the 
exchange of two gluons in the f-channel in the high en- 
ergy limit » 5 where sjs is the square of the four- 
momentum of the colliding hadrons/partons has been 
studied in the Ar^ -factorization approach. In either limit, 
the expansion can be reorganized to simplify the eval- 
uation of the dominant contribution. Away from this 
regime, corrections to each of these approaches may be 
important. Instead of calculating the full NNLO con- 
tributions to the CSM, NNLO* calculations consider 
only real gluon emission at NNLO and control the di- 
vergences using a cutoff. This method has large uncer- 
tainties arising from the sensitivity to the cutoff and the 
renormaUzation scale [18]. We briefly summarize re- 
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cent developments here. For a more complete discus- 
sion, see Refs. [31, 32]. 

2.1. Current progress: Tevatron and RHIC 

The contributions of the NLO color-singlet correc- 
tions reduce the discrepancy between the inclusive CSM 
cross sections and the CDF T data [9]. However, the 
NLO rate still falls too steeply at high pr to describe 
this region successfully. The NNLO* contribution may 
be able to fill the gap between calculations and data at 
high pt- While the large theoretical uncertainties do not 
place strong constraints on the color-octet contribution, 
the data no longer require them [32]. 

Higher-order color singlet contributions to J/ifr and 
i/'(25 ) production have also been calculated. The calcu- 
lation is simpler for the i//(2S) because there is no feed- 
down from higher charmonium states. CDF extracted 
direct if/{2S ) production [9] . The rates were compared to 
the CSM calculations. While the higher order CSM cor- 
rections significantly reduce the discrepancy between 
the calculation and the Tevatron data, they are insuffi- 
cient to remove it entirely. At intermediate pr, the up- 
per limit of the NNLO* calculation is in agreement with 
the data while, for pr > 10 GeV, there is a gap between 
the calculations and the data. The same is true for the 
[32]. 

NLO corrections to '5o and i color octet J/iJ/ [27] 
and T [33] production have been calculated. In both 
cases, these corrections are small for pr < 20 GeV. Val- 
ues of the NRQCD {0^''^Csf'^)} and {0-'^'''C S^^^)} ma- 
trix elements were obtained by fitting the prompt pro- 
duction rate measured by CDF [7] and were found to be 
compatible with those extracted at LO [27]. Feeddown 
was ignored and the P state color octet contribution set 
to zero. Since a satisfactory fit could not be obtained 
for pj < 6 GeV, these points were not included in the 
fit. Note that NRQCD factorization may break down at 
low Pt so that resummation effects may be necessary to 
describe the data in this region. 

Two papers [34, 35] recently appeared with complete 
calculations of NLO a, color-octet production through 
order v'* for the ^^i, '^o, and ^Pj channels. Although 
the two calculations agree at the level of partonic cross 
section, the values of the NRQCD matrix elements ex- 
tracted from the two different fitting procedures are in- 
consistent. In any case, all these computations fail to 
comply with the constraint obtained in [36] by assum- 
ing that the e^e^ — > J/tfr + Xnon-cc rate measured by 
Belle [37] comes only from the color-octet processes: 

< (2.0 + 0.6) X IQ-^ GeV^ 



at NLO in a^. If one keeps in mind that the color-singlet 
contribution -assumed to be zero in the equation above- 
saturates the experimental measurement by Belle, it is 
clear that this upper bound is in fact quite conservative, 
even though this measurement may be affected by the 4 
charged track requirement. The violation of this bound 
by these recent NLO analyses, where color-octets dom- 
inate, should be taken rather seriously, especially be- 
cause e^e^ annihilation may be a cleaner probe of the 
inclusive charmonium production mechanisms. 

On top of its earlier J/iff analyses, PHENIX has 
shown preliminary measurements of the pr dependence 
of the i//(2S) cross section at 200 GeV [14]. This is 
the first measurement of the pr dependence of an ex- 
cited charmonium state at RHIC. PHENIX measured 
the feeddown contribution of the if/{2S ) to the J/ 1// to be 
8.6 + 2.3%, in good agreement with the world average. 

STAR has recently published [13] measurements of 
the J/ip cross section in 200 GeV pp collisions for 
5 < Pt < 13 GeV/c. This greatly extends the pj range 
over which J /iff data are available at RHIC. Although 
PHENIX can trigger at all pr, it has so far been limited 
to Pt below about 9 GeV/c [14] because of its much 
smaller acceptance. 

STAR has also recently presented a measurement of 
the fraction of J/tfr produced from B decays [13]. The 
results appear to be consistent with the CDF measure- 
ment at an order of magnitude higher energy. 

Unfortunately, little is known about J/t// feed down 
from higher charmonium resonances such as the Xc at 
RHIC. These may significantly contribute to the yield 
of y/i/r mesons observed by PHENIX and STAR. For 
instance, it is likely that the fraction of J/iff from Xc 
should be independent of pt, y and V^. This may be 
relevant for the interpretation of the high pT STAR data 
[13] which favor NRQCD over CSM production in cal- 
culations with NRQCD at LO [38] and CSM contribu- 
tions up to NNLO* (see [18]). Both these calculations 
do not include feeddown. At lower pT, PHENIX J/tfr 
data [14, 15] agrees with LO NRQCD including feed- 
down [39]. However, in the latter case, the low pt val- 
ues may call the validity of the perturbative calculation 
into question, especially since the LO NRQCD spec- 
trum diverges as pT 0. Recently, NLO corrections 
to the CSM were evaluated as well and were shown to 
close the gap between the LO CSM and the PHENIX 
data m Pt - I - 2 GeV. When the NNLO* contribu- 
tion is included at pT > 5 GeV, the upper limit agrees 
with data [40]. We note here that the NLO CSM do 
not show any divergences at low pT- However, resum- 
mation of initial state radiation may further improve the 
agreement with PHENIX data. 



4 



2.2. First LHC results 

The first LHC pp run, at 7 TeV, has ah^eady produced 
superior quahty data after a short time. The CMS [41], 
ATLAS [42], ALICE [43] and LHCb [44] experiments 
reported inclusive pr distributions from integrated 
luminosities of 100 nb"', 9.5 nb"', 11.6 nb"' and 14.2 
nb respectively. Before ALICE differential cross sec- 
tions at mid-rapidity are available^, no low pr data 
will be available at mid-rapidity, see Fig. 1. ATLAS 
and CMS experiments are equipped with large mag- 
netic fields, so their single muon pr thresholds are too 
high for reconstruction of low pr J/ip's. However, the 
data taken at pr > 4 GeV/c are consistent with each 
other even though the rapidity ranges differ slightly. 
The CMS rapidity bin, \y\ < 1.4, overlaps the two 
midrapidity ranges reported by ATLAS, |y| < 0.75 and 
0.75 < \y\ < 1.5. 
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Figure 1 : The preliminary mid rapidity inclusive J /if/ pj distributions 
measured by CMS (|y| < 1.4) and ATLAS (|,v| < 0.75 and 0.75 < |y| < 
1.5). Data compilation courtesy of H. Wohri. 

At more forward rapidities, the single muon pj from 
J I ill decay is high enough for reconstruction at 
Pt ^ 0. A compilation of the forward J/tlr data is 
shown in Fig. 2. The CMS and ATLAS data, taken 
in similar rapidity intervals of 1.4 < \y\ < 2.4 and 
1.5 < \y\ < 2.25 respectively, agree rather well for 
Pt > 3 GeV/c. They seem to differ somewhat at lower 
Pt, but the difference is not statistically significant. The 
more forward LHCb and ALICE data agree quite well 



- ALICE has released the inclusive differential pr cross sections 
at large-rapidity, but only the invariant yields are available at mid- 
rapidity at the time of this writeup. 
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Figure 2: The preliminary forward rapidity inclusive J/i// pj distiibu- 
tions measured by CMS (1.4 < ^ < 2.4), ATLAS (1.5 < |y| < 2.25), 
LHCb (2.5 < y < 4.0) and ALICE (2.5 <y < 4.0). Data compilation 
courtesy of E. Scomparin and H. Wohri. 



with the shape of the somewhat more central rapidity 
CMS and ATLAS data. 

The data shown in Figs. 1 and 2 are for inclusive J I if/ 
production which includes feeddown from the Xc and 
ifr(2S) states and semileptonic B meson decays. The 
fraction of J ftp resulting from B meson decays. 



B fraction = 



B J/i/zX 
any inclusive J/ ifr 



(2) 



is compiled in Fig. 3. The preliminary LHC results at 
■\/s - 1 TeV are compared with the high statistics CDF 
midrapidity data [7] at = 1 .96 TeV. The result seems 
to be remarkably independent of both center-of-mass 
energy and rapidity. The lower statistics and lower pj 
STAR data also agree with the magnitude of the higher 
energy data shown in Fig. 3. 

Although the level of the B meson contribution to J/tf) 
production seems to be remarkably independent of V^, 
the average pj of the measured inclusive J/tfr's has been 
seen to increase with V^. The value of both {pj} and 
(Pj) increases approximately linearly with In V* from 
Vs ~ 17 GeV to 200 GeV The preliminary ALICE mea- 
surement reported here [45], seems to follow this trend. 
The dependence of (pj) with is shown in Fig. 4. 
Note that the fixed-target results are centered around 
midrapidity while two rapidity ranges, \y\ < 0.35 and 
1.2 < \y\ < 2.2, measured at RHIC, suggest that the av- 
erage pj is higher at midrapidity than forward rapidity, 
not a surprising result. The ALICE measurement, using 
the muon spectrometer, is also at forward rapidity. 
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Figure 3: The preliminary fraction of non-prompt //lA as a function 
of pt measured by CMS (|.v| < 1.4 & 1.4 < |y| < 2.4), ATLAS {\y\ < 
2.25) and LHCb (2.5 < >' < 4). Data compilation courtesy of H. 
Wohri. 
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Figure 4: The preliminary measurement of (pj.) from ALICE com- 
pared with lower energy data [45] . 



We note that comparisons of these data with model 
calculations using standard DGLAP evolution for the 
gluon distribution with perturbative QCD agree rather 
well (see Fig. 5 for a comparison with the CSM) sug- 
gesting that, even at rather forward rapidity, in a region 
of X and where the parton densities are not well 
measured, no alternative factorization schemes involv- 
ing e.g. saturation effects, are required. 

Finally, CMS also reported a preliminary T(15) mea- 
surement based on 678 + 38 events. The T(15), T(25) 



Figure 5: Comparison between the LO CSM prediction for 
da^jj^" /dy X Br from gg fusion LO contributions in pp collisions 

at VJ = 7 TeV compared to ALICE [46], ATLAS [47], CMS [48] 
and LHCb [49] results multiplied by the direct fraction (and coiTect- 
ing for the non-prompt component, assumed to be 10%, if applicable). 
Adapted from [29]. 



and T(3S ) peaks are clearly separated in the data. 

2.3. What's next ? 

The LHC quarkonium data shown here is just the tip 
of the iceberg. More data will be forthcoming in the 
next runs. Presumably direct J I if/ and T(15' ) production 
will also be measurable, i.e. the feeddown contributions 
will be distinguished and separated. 

On the other hand, after September 2011, the Teva- 
tron will unfortunately be shut down. There are 8.5 fb ' 
of data from Run 11 on tape with ~ 10 fb ' expected 
by the end of the run. This should be enough data to 
measure the Jjij/ and iff(2S) up to pr > 30 GeV/c. The 
three Y S states will also be measured with relatively 
high precision. Not only should the final Tevatron mea- 
surements include the feeddown of Xc to and;^';, to 
T, it is expected that the Xc and Xb P states with large 
branching ratios to the S states can be separated from 
each other using photon conversions. Even though Run 
II will come to an end, the data will continue to be an- 
alyzed for some time with more exciting results forth- 
coming. 

HERA has already ended its physics run. However, 
analyses are still in progress. The HI and ZEUS data 
have been combined in these final analyses, increasing 
the available statistics. These fits will better improve 
our understanding of quarkonium photoproduction. 

Indeed, utilizing the entire Tevatron and HERA data 
sets to study the balance between production of color 
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singlet and color octet states in hadro- and photopro- 
duction will provide more insight into the production 
mechanisms than fitting to either data set alone. Last but 
not least, the B-factory analyses on charmonium pro- 
duction, not only from B decays, are detailed enough 
nowadays to separate out the different inclusive contri- 
butions, imposing constraints on the production mecha- 
nisms at B factories, ep colliders and pp colliders. For 
instance, it starts to be rather clear that C = +1 color- 
octet transitions are much less probable than initially 
thought, which suppresses the leading color-octet yield 
at pt < 5 GeV in hadroproduction. 

In 200 GeV pp collisions, PHENIX has measured 
cross sections for production of the (unresolved) T 
states and the i/'(25). A Xc measurement will be pub- 
lished soon. All of these measurements involve low 
yields and will benefit greatly from improved luminos- 
ity in the next few years. PHENIX has also measured 
the polarization of the in 200 GeV pp collisions and 
will also do so with existing 500 GeV data. PHENIX 
and STAR will also produce unpolarized J/il/ and T re- 
sults at 500 GeV, a new energy midway between the 200 
GeV data and the Tevatron energy. 

RHIC is vigorously pursing a long term plan of de- 
tector and machine upgrades. The ongoing RHIC lumi- 
nosity upgrades will be completed by the 2013 run. The 
introduction in 2011 and 2012 of silicon vertex detec- 
tors into PHENIX at mid and forward rapidity and of the 
STAR Heavy Flavor Tracker in 2014, will enable open 
charm and open bottom to be measured independently 
with greatly improved precision. The detector upgrades 
will also allow improved measurements of quarkonium 
states due to improved mass resolution and background 
rejection. 

The RHIC run plan for the next 5 years or so will 
be centered on exploiting the capabilities of the new 
silicon vertex detectors and other upgrades, combined 
with the increased RHIC luminosity. There will likely 
be long pp runs at 200 GeV, plus shorter runs at 62 
GeV to explore the energy dependence of open and hid- 
den heavy flavor production. The luminosity increase 
and the PHENIX detector upgrades will enhance the 
PHENIX heavy quarkonium program with increased pr 
reach for the J/if/ and low statistics measurements of the 
combined T states. The increased luminosity will en- 
able the large acceptance STAR detector to extend its 
Pj reach to considerably higher values for T and J/iJ/ 
measurements than previously possible. 

The RHIC program for the period beyond about 5 
years is still under development. The RHIC experiments 
are presently engaged in preparing a decadal plan that 
will lay out their proposed science goals and detector 



upgrades for 201 1 to 2020. 

3. Quarkonium production in pA collisions 

The study of charmonium production is an interesting 
test of our understanding of strong interaction physics. 
For pp collisions, various theoretical approaches have 
been proposed, but a satisfactory description of J/tfr 
hadroproduction is still missing [31]. In spite of this 
situation, the study of quarkonium production in the ap- 
parently more difficult proton-nucleus environment has 
also attracted considerable interest. In fact, in pA, the 
heavy-quark pair is created in the nuclear medium, and 
the study of its evolution towards a bound state can add 
significant constraints to the models. For example, the 
strength of the interaction of the evolving cc pair with 
the target nucleons, that can lead to a break-up of the 
pair and consequently to a suppression of the J/t// yield, 
may depend on its quantum states at the production level 
(color-octet or color-singlet), and on the kinematic vari- 
ables of the pair [50, 51]. In addition to final-state 
eff'ects, also initial-state effects may influence the ob- 
served J/ip yield in pA. In particular, parton shadowing 
in the target nucleus [52, 53, 54, 55, 56] may suppress 
(or enhance, in case of antishadowing) the probability of 
producing a J/ 1/', and its effect depends on the kinemat- 
ics of the cc pair production [57]. Moreover, the energy 
loss in the nuclear medium of the incident parton [58], 
prior to cc production, may significantly alter the J/ ij/ 
cross section and its kinematic distribution. Finally, a 
suppression of the Jl^ has been proposed a long time 
ago as a signature of the formation, in ultrarelativistic 
nucleus-nucleus collisions, of a state where quarks and 
gluons are deconfined (Quark-Gluon Plasma) [59]. Re- 
sults from pA collisions, taken in the same kinematical 
conditions of AA, and properly extrapolated to nucleus- 
nucleus collisions, can be helpful to calibrate the contri- 
bution of the various cold nuclear matter effects to the 
overall observed suppression [60, 61]. 

3.1. Basic phenomena and observables 

3.1.1. Hierarchy of titne scales 

Nuclear effects are controlled by the characteristic 
time scales of charmonium production. A colorless cc 
pair created in a hard reaction is not yet a physical char- 
monium state since it does not have a fixed mass. Ac- 
cording to the uncertainty relation it takes time to disen- 
tangle the ground state charmonium from its excitations. 
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This time scale can be also interpreted as formation time 
of the charmonium wave function. At charmonium en- 
ergies Ecc < 25 GeV in the nuclear rest frame this time 
is shorter than the mean nucleon spacing in a nucleus, 
tf < 2fm, so one can treat the formation process as 
instantaneous, and the attenuation in nuclear matter is 
controlled by the inelastic charmonium-nucleon cross 
section, 

5(L) = exp[-(rf^^Lp^], (4) 

where L is the path length in the medium of density p^- 
At much higher energies, Ecc » 100 GeV, the forma- 
tion time is long, // » Ra, even for heavy nuclei, the 
initial size of the cc dipole is "frozen" by Lorentz time 
dilation during propagation through the nucleus, so the 
attenuation factor gets the simple form. 



S{L) = <exp [-o-{rT)LpA\),^ 



(5) 



where the exponential is averaged over transverse sepa- 
ration rj of the dipole weighted with the initial and fi- 
nal distribution amplitudes [50]. This regime is relevant 
for JI4i production at RHIC at the mid rapidity, where 
Ecc - 300 GeV. The energy rises with y as e' rela- 
tive to the target nucleus, but decreases as e"^ relative to 
the beam nucleus. Therefore, the asymptotic regime is 
relevant to any positive rapidity in dA collisions, but in 
AA collisions at forward rapidities nuclear effects turn 
out to be in different regimes for the two colliding nu- 
clei. 

Another important time scale characterizes the hard 
process of cc pair creation. Although the proper time is 
short, f ~ l/rric, it is subject to Lorentz time dilation in 
a frame, where the charm quarks have high energy. 



2Ec, 
Ml 



(6) 



If this time is short, fe < 2 fm, one can consider pro- 
duction as instantaneous. This regime corresponds to 
charmonium energy Ecc < 50 GeV. This estimate is rel- 
evant for production of x, but the characteristic energy 
may be about twice higher for J/t//, because in the color 
singlet model the invariant mass of cc is (M?^) = '^■iii]/^- 
This time scale is usually short for J/tfr produced at 
mid rapidity in fixed-target experiments. For instance 
at = 40 GeV tc = 1.2 fm if = 0. 

The most difficult for calculations is the intermediate 
energy regime, where either tc or/and tf are of the order 
of the nuclear size. Theoretical tools for this case, the 
path-integral technique, was developed in [50, 62]. The 
cc dipole is produced with a separation ~ l/nic and ends 
up with the large size of J /iff. The time scale for such 



an expansion, tf, rises with energy, therefore the effec- 
tive absorption, or break-up, cross section is expected to 
drop with energy [50]. Indeed, this effect was observed 
in fixed target experiments, as is plotted in Fig. 6. Data 
are well explained by color transparency for expanding 
cc dipoles. 
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Figure 6: The effective break-up cross section as function of // {fj en- 
ergy. The solid curve shows theoretical expectation [62]. The dashed 
curves show the estimated theoretical uncertainty. 



3.1.2. Higher twist charm quark shadowing 

In the regime of long time scale for cc production, of 
tc » Ra, the process, g — > cc, is subject to shadow- 
ing. This is apparently higher twist, since the mean cc 
transverse separation is {rj) ~ l/nic. Same is true for 
the break-up cross section of the produced colorless cc 
propagating through the nucleus. 

The amplitude of cc production at a point with im- 
pact parameter b and longitudinal coordinate z, aver- 
aged over the dipole size, reads [63, 64], 



SpA(b, z) 
1 



X exp 



d^rr W-ccirr) 



1 



(7) 



:0--ccArT)T-{b,z) - -o--ccirT)T+{b,z) 



}iemTAb,z) = £^dz'pA(b,z'y, TAb,z) = T^W - 
T-(b,z), and TAib) - T-{b,oa). Shadowing for cc pro- 
duction over the nuclear thickness T-(b,z) occurs with 
the shadowing cross section corresponding to a 3-body 
dipole, gluon and cc, which for equal momenta of c and 
c equals to o-ccgirr) = lo-cc(rT/2) - lo-cdrr)- 
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For tc » Ra the weight factor in (7) has the form 
[63], Wccirr) ^ Koimcrr) 4 i'JI^rr), where one fac- 
tor kt comes from the amplitude of cc production, and 
another one either from the amplitude of gluon radia- 
tion in the case production, or from the radial wave 
function of xi- 

With the survival probability amplitude Eq. (7) the 
nuclear ratio reads. 



oo 



(8) 



The results at = 200 GeV are plotted by dashed 
curve in Fig. 7 as function of y. [64]. 




Figure 7: Dashed curve presents nuclear suppression of J/i// as func- 
tion of rapidity in pA collisions. Solid curve is corrected for gluon 
shadowing. Data are for dAu collisions at V* = 200 GeV [65]. 



hadronic representation [68, 69] and in the global anal- 
yses EKS98 [53, 54], EPS08 [55] and EPS09 [56]. We 
emphasize that the two last analyses included data on 
proton-nucleus collisions, which makes them model de- 
pendent. We finally note that nDSg and EKS98 param- 
eterizations are compatible for x < 10"^. 

Also, the choice of the adequate partonic production 
mechanism - either via a2^1ora2^2 process - 
may affect both the way to compute the nuclear shad- 
owing and its expected impact on the production [57]. 

This controversy should be settled after nuclear ef- 
fects for charmonium production at LHC will be studied 
experimentally. 

3.1.4. Non trivial transition from pA to AA 

The usual strategy in search for a signal of final state 
attenuation of charmonia in the created dense medium, 
is extrapolation to AA collisions of the cold nuclear mat- 
ter effects observed in pA. This is how the "anomalous" 
suppression of was observed in the NA50/60 data 
[70, 71, 72, 73, 74]: one fits the effective absorption 
cross section cToii to describe the nuclear suppression 
observed in pA collisions. Some results of such a pro- 
cedure are shown in Fig. 6. Then with this cross section 
one predicts the normal cold nuclear matter effects in 
nuclear collisions, and whatever deviates from such an 
expectation is called anomalous and is related to the fi- 
nal state suppression. 

However, recent data [75] from the NA60 experiment 
for pt broadening of 7/ 1/' depicted in Fig. 8 demonstrate 
that such an extrapolation may not be straightforward, 
since the "cold" nuclear matter in AA collisions turns 
out not to be cold at all. 



3.1.3. Leading twist gluon shadowing 

In terms of the Fock state decomposition gluon shad- 
owing corresponds to higher Fock components in the 
projectile hadron, containing at least one gluon. Mul- 
tiple interactions of this gluon in the nuclear target are 
the source of shadowing. 

Theoretical expectations for gluon shadowing are 
quite diverse. A weak shadowing was predicted within 
the dipole approach [66] and in an analysis [52] of DIS 
data developed at leading and next-to-leading order - 
nDS LO and NLO parameterizations^^ -. However, a 
rather strong shadowing and correspondingly antishad- 
owing for gluons was evaluated in [67] based on the 
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Figure 8: J/if/ pr broadening measured at SPS as a function of the 
path length L in nuclear matter 



'We note that the same DIS analysis provided a constrained NLO 
fit with a stronger gluon shadowing - nDSg -. The of this fit is 
however not as good as that of the unconstrained fit -nDS. 



Indeed, broadening, Ap^, is known to be proportional 
to the path length in the medium and its density. The 
data show that broadening observed in nucleus-nucleus 
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Figure 9: Left: A two-dimensional (time - longitudinal coordinate) plot for charmonium production in a collision of nuclei A and B in the cm. 
of the colliding nucleons. The solid and dashed lines show the nucleon and gluon trajectories, respectively. Right: Nuclear suppression of Jjifj 
production in pA and AB collisions as function of the product Ax B. The two curves demonstrate theoretical uncertainty of calculations [76] . The 
circles and squares correspond to pA and AB data, respectively. The data points are from [70,71,72]. 



collisions is twice as large as in pA for the same total 
path length. This means that the cold nuclear matter is 
twice as dense in AA compared to pA collisions. Cor- 
respondingly, the effective absorption cross section a-ghs 
should be taken twice bigger when one makes predic- 
tions for AA collisions. Such a significant enhancement 
of absorption effects may explain the observed anoma- 
lous suppression. 

This effect was successfully predicted in [76] as a re- 
sult of multiple A^A^ collisions and gluon radiation pre- 
ceding their interaction with the Jlij/ (or cc dipole), as is 
illustrated in the left panel of Fig. 9. 

It was also shown that quantitatively this effect is able 
to explain the observed anomalous J I if/ suppression, as 
is depicted in the right panel of Fig. 9. Simultaneously, 
the same effect of increased medium density leads to a 
significant increase of broadening in AA compared to 
pA collisions. 

Suppression of Jj^ caused by gluon radiation de- 
creases as 1/ ^fs and vanishes at high energies of RHIC 
and LHC. However, at the energies of LHC broadening 
of J/i// is expected to be highly enhanced in AA com- 
pared to pA collisions for another reason. 

Notice that observation of an anomalously large 
broadening has an important advantage compared to 
anomalous suppression. While the latter may originate 
from either initial, or final state interactions, which can 
be easily mixed up, the former comes entirely from ini- 
tial state multiple interactions. Indeed, a colorless cc 
dipole attenuates in a medium, but does not have any 
energy loss and does not change its momentum. Note 
that the partonic process may affect this result [77]. 



3.1.5. Transverse momentum broadening 

The mean transverse momentum squared of heavy 
quarkonia increases in pA collisions compared io pp. 
The magnitude of the effect, A;?^ - {p\)a - {Pj}p is 
usually called broadening. Broadening of a parton prop- 
agating through a nuclear medium can be calculated 
within the dipole approach as [78, 79], 



2C{E) dzpAiz), 



(9) 



where the factor C(E) is related to the dipole-proton 
cross section cr^^(rj) known from phenomenology. For 
broadening of gluons 



rT=0 



(10) 



Using the dipole cross section in the saturated form fit- 
ted to DIS data [80, 66] one can predict the factor Cg(E) 
and broadening Eq. (9). The results depicted in Fig. 10 
by dashed curves well agree with data for J/i// and T. 
Inclusion of gluon shadowing corrections reduces the 
amount of shadowing, but within the error bars. 

One should also understand the pj dependence of the 
cross section and the energy dependence of the mean 
value (pj), depicted in Fig. 11. 

Remarkably, data on pp, pA and even AA collisions, 
at the energies of fixed target experiments and at RHIC 
[65], are described well by the simple parametrization. 



dcr 
dpj 



1 + 



6{pl}) 



(11) 



10 




2 Figure 12: The exponent a = I + ln(i?p^)/ ln(A) as function of pr 

calculated with Eq. (1 1) in comparison with data from [82]. 



Figure 10: Broadening for J/i// and T [81] is shown by circles and 
triangles respectively. The dashed and solid curves coiTespond to the 
predictions based on Eq. (9) without and with the corrections for gluon 
shadowing respectively. 
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Figure 11: The J/ip {p^) in pp collisions as function of energy. Data 
are from [10]. 



3.1.6. Feynman xp dependence of nuclear effects 

Available data on jc/r-dependence of nuclear effects in 
charmonium production demonstrate increasing nuclear 
suppression towards large Xf, as is shown in Fig. 13, for 
the exponent a{xf) characterizing the A-dependence, 
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If one considers that broadening does not alter the shape 
of the pr-distribution of produced J/tfr, the simplest way 
to calculate the pj-dependence of the pA cross section 
is just making a shift Apj in the mean value {pj} for pA 
compared to pp. The resulting A-dependence of the pA 
over pp ratio, RpA, is presented in Fig. 12 as function 
of pt in comparison with data [82] at = 39 GeV. 
The A-dependence is parametrized as RpA - A'' and 
a(pj) is plotted in Fig. 12. 



Figure 13: The exponent a(Xf) characterizing the A-dependence, 
^" at = 39 GeV [82] an = 20 GeV [83]. 

Although several mechanisms contributing to the 
J/ifr suppression have been proposed, there is still no 
comprehensive description so far of the observed xp- 
dependence of the nuclear effects. The first mechanism 
proposed in [84] explored initial state energy loss in- 
dependent off the collision energy, according to pQCD 
calculations, or string model. However, a finite energy 
loss becomes insignificant at high energies, while data 
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in Fig. 13 demonstrate an approximate Xf scaling. Ap- 
parently, the energy loss must be proportional to the ini- 
tial energy, in order to provide an x/r-scaling. Such an 
energy loss was ad hoc introduced in [58], and was crit- 
icized later in [85], as contradicting perturbative QCD 
calculations. 

The observed -dependence of nuclear suppression 
was interpreted in [86] in terms of the Fock state decom- 
position as the energy sharing problem enhanced by the 
nucleus towards the kinematic limit. This mechanism 
also explains a similar suppression observed at large xp 
in many other reactions, Drell-Yan process, high and 
low pt production of different species of hadrons, etc. 

Attempts to explain the observed suppression to glu- 
ons shadowing, as a dominant mechanism, was not suc- 
cessful. That would lead to xi scaling, which is severely 
broken in data [81]. 

3.2. Existing measurements and future perspectives 

3.2.1. Fixed-target experiments 

Since in charmonium production in pA a rather com- 
plicated interplay of various physical processes occurs, 
the availability of accurate sets of data, spanning large 
intervals in the incident proton energy, and covering 
large xp and pj regions, is essential for a thorough un- 
derstanding of the involved mechanisms. At fixed tar- 
get energies, high-statistics Jlij/ samples have been col- 
lected in recent years by the DESY experiment HERA- 
B [87], at 920 GeV incident energy, by E866 [88] at 
FNAL at 800 GeV, by the CERN-SPS experiment NA50 
at 400 and 450 GeV [89] and by the NA60 experi- 
ment at 158 and 400 GeV [90]. Usually, nuclear effects 
have been parametrized by fitting the A-dependence of 
the production cross section with the simple power law 



^'Tjifi ~ "^y/iA ' ^^'^ studying the evolution of 

a with Xp and p-p. Alternatively, nuclear effects have 
been expressed by fitting the data in the framework of 
the Glauber model [91], having as input parameters the 
inelastic nucleon-nucleon cross section and the density 
distributions for the various nuclei [92]. The model 
gives as an output the so-called absorption cross 
section cr^^^. Clearly, both a and cr^*^ represent effec- 
tive quantities, including the contribution of the various 
sources of nuclear effects sketched in the previous sec- 
tion. 

The results for a as a function of xp are summarized 
in Fig. 14, and show various remarkable features. First 
of all one can note a steady increase in the strength of 
nuclear effects {a decreases) when increasing xp. More 
in detail, the decrease of a at large xp has usually been 
attributed to energy loss of the projectile partons and 



indeed various models have tried to explain, at least 
qualitatively, such an effect (see e.g. [93] and more re- 
cently [77]). On the other hand an important constraint 
on the size of parton energy loss comes from experi- 
ments studying nuclear effects on Drell-Yan production. 
At this workshop [94], it has been shown that the initial- 
state parton energy loss needed to describe the NA3 data 
on Drell-Yan production [95] is too small (by a factor 5 
to 10) to explain at the same time the observed high- 
xp suppression for Therefore, a quantitative as- 

sessment of this effect is still lacking. At the other ex- 
treme of the explored xp region (negative values), one 
observes a slight tendency towards enhancement of the 
y/i/' production (a> 1). In this region a model [96] for- 
mulated in the framework of reggeon phenomenology, 
indeed predicts an anti-screening effect, which cannot 
be reproduced by more conventional approaches based 
on the commonly used shadowing parameterizations. 
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Figure 14: Compilation of experimental results on vs Xf . 

At midrapidity, several experiments have studied Jj\p 
production. The most remarkable feature is a tendency 
towards stronger nuclear effects when of the pA 
interaction decreases. In this region, the relevant pro- 
cesses affecting the production should be the final-state 
breakup of the cc pair and the parton shadowing in the 
nuclear target. Both processes should approximately 
scale with X2, i.e. one would expect nuclear effects to 
be the same at a certain X2 independently of the initial 
proton energy. However, recent results by NA60 [97], 
comparing or as a function of xo for y/t/r production at 
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158 and 400 GeV, show that such a scaling does not 
hold, and that therefore other mechanisms should be in- 
voked to explain the V^-dependence of nuclear effects 
at midrapidity. 
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number of A^A^ collisions, is formed. There is a remain- 
ing systematic uncertainty due to the centrality depen- 
dence of the tracking and particle identification efficien- 
cies. 
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Figure 15: Aq;=Q'4oogcV ~ flssGcV vs .V2, as measured by NA60 in pA 
collisions at 158 and 400 GeV. 

To summarize, in spite of the rather extended set of 
measurements, which clearly define a trend for nuclear 
effects on J/iJ/ production, the interpretation of fixed- 
target observations still remains unsatisfactory. For the 
future, further results obtained at lower energies, at the 
SPS or at the future FAIR facility, may help produce 
a better characterization of the evolution of nuclear ef- 
fects vs ^fs. First ideas for fixed-target measurements at 
higher yfs [98], using proton beams extracted from the 
LHC, are also under consideration. 

3.2.2. Collider experiments 

The study of J/ tfr production and propagation in cold 
nuclear matter has also been carried out at RHIC en- 
ergies, through dAu collisions. The conclusions from 
the first round of measurements from the PHENIX ex- 
periment were not so sharp [65], because of the limited 
available statistics, but this problem has been recently 
overcome, with the availability of a first set of results 
corresponding to a much larger data sample (Run-6). 
One of the key features of the PHENIX results is the 
large rapidity coverage, including positive and negative 
rapidities, corresponding to a large range of relative mo- 
menta between the cc pair and the medium. 

More in detail, PHENIX has recently released c/Au 
Rcp data for J/t// production [14] in nine rapidity bins 
over \y\ < 2.4. Systematic uncertainties associated with 
the beam luminosity, detector acceptance, trigger effi- 
ciency, and tracking efficiency cancel when Rcp, the 
ratio of central to peripheral events, normalized to the 
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Figure 16: The elfective absorption cross section as a function 
of rapidity extracted from PHENIX dAu Rcp data in the extrinsic 
2 — * 2 scheme [102] (in red closed circle) compared to the 2 — > 1 
scheme [99] (in blue open circle) using a) EKS98, b) EPS08 and c) 
nDSg. From [102]. 



More importantly, there are significant systematic un- 
certainties in the centrality dependence of Rcp due to the 



13 



use of a Glauber model to calculate the average number 
of nucleon-nucleon collisions as a means of estimating 
the relative normalization between different centrality 
bins. The systematic uncertainty due to the Glauber cal- 
culation is independent of rapidity. 

The PHENIX dAu Rqp data have been independently 
fitted at each of the nine rapidities [99] employing a 
model including shadowing and Jjij/ absorption. The 
model calculations [100] use the EKS98 and nDSg 
shadowing parameterizations with < cr^*^ < 15 mb. 
The best fit absorption cross section was determined at 
each rapidity, along with the +lcr uncertainties asso- 
ciated with a) rapidity-dependent systematic uncertain- 
ties and b) rapidity-independent systematic uncertain- 
ties. The most notable feature is the stronger effective 
absorption cross section at forward rapidity, similar to 
the behavior observed at lower energies [88]. In fact, it 
is striking that the extracted cross sections at forward 
rapidity are very similar for PHENIX ( ^[s^ - 200 
GeV) and for the E866 collaboration ( ■yjs;;'^ = 38.8 
GeV) [101], despite the large difference in center-of- 
mass energy. One should note the large global system- 
atic uncertainty in (Tabs extracted from the PHENIX Rcp 
data, dominated by the uncertainty in the Glauber esti- 
mate of the average number of collisions at each central- 
ity. Although it does not affect the shape of the rapid- 
ity dependence of o"}^^, it results in considerable uncer- 
tainty in the magnitude of the effective absorption cross 
section. 

It was recently suggested [57, 102] that the large in- 
crease in effective absorption cross section at forward 
rapidity obtained from a LO CEM calculation [99], may 
be significantly moderated if a more accurate 2 — > 2 
kinematics is used. Figure 16 shows the results from 
the 2 — » 1 kinematics [99], compared to the ones from a 
2^2 kinematics [102] using the i-channel cut [103] as 
partonic model. The difference emphasizes the impor- 
tance of understanding the underlying production mech- 
anisms. 

PHENIX has very recently released [104] final R^au 
and Rcp data from the 2008 liAu RHIC run (where RdMi 
is the y/i/r yield in dAu collisions normalized to the pp 
yield times the number of A^A^ collisions). The final Rcp 
data are in very good agreement with the preliminary 
results described above. A comparison of the RdAu data 
with the Rqp data shows that a simultaneous description 
of the two observables will require a stronger than linear 
dependence of the Jjijj suppression on the nuclear thick- 
ness function at forward rapidity. The dependence of the 
suppression on nuclear thickness is at least quadratic, 
and is likely higher This result may have important 
implications for the understanding of forward-rapidity 



dPM physics, as well as for the estimate of cold nuclear 
matter effects in AuAu collisions. 

3.2.3. Cold nuclear matter effects at LHC energy 

The new energy regime recently opened up by the ad- 
vent of the LHC offers new possibilities for the study 
of quarkonium in cold nuclear matter. However, al- 
though the study of pA collisions is foreseen in the LHC 
physics program, the details of the acceleration scheme 
and the choice of the nuclear beam still have to be final- 
ized. Consequently, the study of the physics observables 
accessible to the experiments is still at a rather prelim- 
inary level. Clearly, pA data at the LHC are essential 
to investigate various physics effects, and "in primis" 
nuclear shadowing in a still unexplored x-range. The 
final-state interaction of charmonia with cold nuclear 
matter is another attractive topic. Due to the expected 
very short overlap time of the heavy-quark pair with nu- 
clear matter, one could expect its influence on the still 
almost point-like cc to be very small. However, these 
considerations are very qualitative and deeper theoreti- 
cal studies are clearly necessary. 

At this workshop, preliminary studies carried out in 
the frame of the ALICE experiment have been dis- 
cussed [105]. In particular, the availability of a forward 
muon spectrometer (rapidity coverage 2.5 < y < 4) 
gives the opportunity for studying effects connected 
with gluon saturation via the measurement of Rpp\, for 
the J/t//. A model for heavy-quark production in a 
Color-Glass Condensate environment [106] has been 
used to make predictions for Rpp\, as a function of y and 
Pj in the acceptance of the ALICE muon spectrometer. 
Such predictions quantitatively differ from the expec- 
tations related to a pure shadowing scenario, making 
such a measurement an interesting tool to investigate 
saturation-related effects. 

4. Polarisation and feed-down in quarkonium pro- 
duction 

Until recently, the numerous puzzles in the predic- 
tion of quarkonium-production rates at hadron collid- 
ers were attributed to non-perturbative effects associ- 
ated with channels in which the heavy quark and an- 
tiquark are produced in a colour-octet (CO) state [107, 
108, 109, 110]. It is now widely accepted that and 
al corrections to the CSM [3] are fundamental for un- 
derstanding the pt spectrum of J/t// and T produced in 
high-energy hadron collisions [24, 23, 25, 26, 18, 32]. 
The effect of QCD corrections is also manifest in the 
polarisation predictions. While the J/i/f and T produced 
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inclusively or in association with a photon are predicted 
to be transversally polarised at LO via a color singlet 
(CS), it has been recently emphasised that their polari- 
sation at NLO is increasingly longitudinal when pj gets 
larger [25, 18, 127, 128]. 

On the other hand, the LO NRQCD calculation (for 
which CO transitions dominante) predicts a sizable 
transverse polarisation rate for large pj J/t// [111, 112, 
1 13, 1 14, 115] whereas the Tevatron CDF measurement 
at Fermilab [8] displays a slight longitudinal polarisa- 
tion at large pr for J/tfr and Y and a stronger one for 

To clarify the situation, more efforts on both exper- 
imental and theoretical aspects are expected from the 
forthcoming LHC. In the following, we review impor- 
tant aspects of quarkonium polarisation both from a 
theoretical and experimental perspective with particu- 
lar emphasis on the J - I states such as the if/{2S) 
and T{nS) . 

4.1. Frames and polarizations 

The polarisation of the J - \ quarkonia is defined 
from their dilepton decay. In general, the angular distri- 
bution of the dilepton is 



1 



(3 + A,,) 



{\+A§co5 i^-H/l^sin i?cos2i^ 



+A§^ sin 2§ cos ^fl), 



(12) 



where § and ip are the (polar and azimuthal) angles of 
the positive lepton with, respectively, the polarisation 
axis z and the production plane xz (containing the col- 
liding particles and the decaying meson). Several polar- 
isation frames have been used in the past. In the helicity 
frame the polar axis coincides with the flight direction 
of the meson in the centre-of-mass frame of the collid- 
ing hadrons. In the Collins-Soper frame, the polar axis 
reflects, on average, the direction of the relative velocity 
of the colliding partons. 

To clarify this puzzling situation, improved measure- 
ments are needed. So far, most experiments have pre- 
sented results based on a fraction of the physical in- 
formation derivable from the data: only one polarisa- 
tion frame is used and only the polar projection of the 
decay angular distribution is studied. These incom- 
plete results prevent model-independent physical con- 
clusions. Moreover, such partial descriptions of the ob- 
served physical processes reduce the chances of detect- 
ing possible biases induced by insufficiently controlled 
systematic effects. In the forthcoming LHC measure- 
ments, it is important to approach the polarisation mea- 



surement as a multidimensional problem, determining 
the full angular distribution in more than one frame. 

In a series of recent works [116, 117, 118, 119], a 
frame-invariant formalism was proposed to minimize 
the dependence of the measured result on the experi- 
mental acceptance. It provides a much better control 
of the systematic effects due to detector limitations and 
analysis biases. The overall idea is to study both exper- 
imentally and theoretically the quantity 



A^ + 3/1^ 
1 - ^ 



(13) 



instead of and /l^. 



4.2. Charmonium polarization 

The NLO QCD corrections to y/f/r polarisation in the 
CSM at Tevatron and LHC have been calculated for 
the first time in [25]. At order a^, one should also 
take into account the contribution for gg — > Jlif/ + 
cc [23]. These are dominant at large pr- The results 
showed that the polarisation changes drastically 
from strongly transverse at LO into strongly longitu- 
dinal at NLO. While, this does not completely agree 
with the CDF data, the trend is rather encouraging es- 
pecially since feed-down from P-waves can alter the 
prompt yield polarisation. The same trend is observed 
for the ^{2S ) [32], for which polarization measurements 
are not affected by the P-waves. This is illustrated by 
Fig. 17. 
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Figure 17: Polarisation of tjj{2S ) directly produced as function of its 
transverse momentum pj- at the Tevatron as predicted by the CSM at 
LO, NLO^ (including gg Jjifjcc) and at NNLO* compared to the 
CDF measurement at V* = 1-96 TeV. 

The first NLO QCD corrections of the Jjxp produc- 
tion via CO states '^S^^^\ at the Tevatron and LHC 
were studied in [27]. Contrary to the CS case, these 
corrections only slightly change the pj distribu- 
tions and the polarisation compared to the LO NRQCD 
calculations. This result implies that the perturbative 
QCD expansion quickly converges for production 
via the S -wave CO state, in contrast to the CSM, where 
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Figure 18: Direct-7/i^ yield polarisation a based on the work [27] 
and the NRQCD matrix elements are from the recent work [34]. The 
experimental data is from the CDF measurement at the Tevatron [8] 
for the prompt J/ip. 

the NLO contributions open new channel which show a 
leading pr behaviour with a different J ftp polarisation. 

As shown in Fig. 4.2, an obvious gap between the 
theoretical results for J/ip polarisation calculated up to 
NLO including CO transition [27] and the experimental 
measurements at Tevatron is observed for increasing pr- 
There still remains only a narrow window which might 
fill the gap if one follows the a priori rigourous theo- 
retical framework of NRQCD, namely the NLO correc- 
tions to J/ 1// production via f -wave CO states and J/ ij/ 
production by feed down from Xcj- It is worth noting 
however that the latter possible solution does not apply 
for 1^(25 ) where a similar gap is observed -with a large 
magnitude even-, while the former one would be at odds 
with the suppression of C = + 1 CO transition expected 
from e^e" studies. 

The results given by two recent works [34, 35, 120] 
show that the pr distribution of the P-wave CO contri- 
bution at QCD NLO can be decomposed into the linear 
combination of the two 5 -wave COs. In Ref. [34], the 
Xc feed down to prompt J/tfr hadroproduction is taken 
into account, contributing at the level of 20 - 40%. In 
this case, the NRQCD matrix elements {O^} are deter- 
mined to be {Ol(^S i)> = 0.0005 GeV^ and (OlCSo)} = 
0.076 GeV^'. It should be noted that some experimental 
data have been omitted since the fit of the pj- distribu- 
tion starts at 7 GeV. The theoretical prediction for the 
polarisation of the direct y/f/r hadroproduction is plotted 
in Fig. 18, based on the calculation of Ref. [27]. We also 
note that such a value for {O'^i^So)} does comply with 
the constraint of Eq. (1). This either hints at an overesti- 
mation of this transition in these studies or at a breaking 
of NRQCD matrix-element universality for charmonia. 
In such a case, NRQCD would become nearly unpredic- 
tive as far as charmonium production is concerned. 

Two aspects should be emphasized when comparing 
the theoretical prediction and the existing experimen- 
tal measurements. It is impossible to give the polarisa- 



tion of the prompt J/if/ hadroproduction without a cal- 
culation of the polarisation for the Xc feed-down con- 
tribution if it really contributes 20 - 40% to the pr 
distribution. In addition, the polarisation prediction in 
Fig. 18 for direct J/i// production assumes that the CO 
^P^j\j/il/) NRQCD matrix element can be neglected. If 
not, its contribution to J/ip polarisation must calculated 
at NLO accuracy. At LO, the Xc ™d the P-wave CO 
contributions were both taken into account in the first 
NRQCD-based predictions [113]. 

4.3. Bottomonium polarization 

As regards the polarisation of Y produced pp col- 
lisions, the NLO CSM correction was first presented 
in [26], and followed by an analysis including the real- 
emissions at NNLO [18] which we discuss later on. The 
polarization is quite similar to that of J I if/, i.e. the polar- 
isation changes drastically from strongly transverse at 
LO into strongly longitudinal at NLO. The NLO results 
are also rather similar to those obtained in the kj fac- 
torisation approach at LO [121, 122], 
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Figure 19: Transverse momentum distribution of polarization param- 
eter a for direct T production at the Tevatron. 

The NLO QCD corrections of the Y production via 
CO [ ' 5 [j*^ V5 f ] at the Tevatron and the LHC were stud- 
ied for the first time in [33]. The Y polarization param- 
eter, a, from CO ***^] as function of pr is shown in 
Fig. 19. Slight changes can be observed when the NLO 
corrections are taken into account. For any pj, the di- 
rect Y yield from CO is transversally polarized, with 
a > 0.5 as soon as pr > 10 GeV. The predictions for 
the polarization of Y production when the CSM NLO 
yield is added are also presented in the figure as a "to- 
tal" result. 

In Fig. 20, the polarization of inclusive (i.e. prompt) 
Y production at the Tevatron is shown. As the polariza- 
tion of Y from the feed-down of XbinP) is not available 
yet, a huge band is obtained only by imposing that the 
polarization of this part is between -1 to 1. The experi- 
mental data from DO and CDF is also shown in the same 
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Figure 20: Extrapolated prompt T(15) polarization as a function of its transverse moinentum at the Tevatron as predicted by NRQCD (including 
CO transitions) [33], CSM at NLO [18] and CSM at NNLO* [18], compared to the DO [22] and CDF [123] data. 



figure. Since both data sets clearly disagree, it is diffi- 
cult to draw any strong conclusion from such a compar- 
ison. Yet, if one follows the CDF analysis based on a 
higher integrated luminosity and where the 3 T states 
are clearly separated out, one observes the same gap 
opening with pr as for the i^'s, even with such a large 
theoretical uncertainty. For the sake of completeness, 
one should underline that the direct yield from CO tran- 
sitions is always predicted to be transverse. The CO 
yield can only become slightly longitudinal thanks to 
the (unknown) feed-down. This is a clear motivation 
for the CDF and DO to measure the polarisation of di- 
rect Y(15) or that of T{3S), the latter being de facto 
direct. 

In addition to the evaluation of the NLO CSM con- 
tributions to T production, a subset of the NNLO con- 
tributions were evaluated in [18], these are the real- 
emission contribution at order believed to provide 
the leading contribution at large pj. It was indeed 
found that these contributions, included in the NNLO* 
yield, could be large and bring an agreement between 
a CSM-based evaluation and the CDF data. It also oc- 
cured that the NNLO* yield showed the same polarisa- 
tion than the NLO yield, excepted for the -important- 
fact that the yield was correctly described. The polar- 
ization of directly-produced Y predicted likewise does 
not agree with the experimental data since it is strongly 
longitudinally polarized. Nevertheless, as we have 
discussed above, the experimental data includes con- 



tributions from P-wave decays, up to roughly 40 %, 
and these can decay into transversally polarised Y(15). 
Even if the feed-down yield is unpolarized, the prompt 
polarization would only be 0.6 times the direct polar- 
ization. Clearly when this is taken into account, the 
CSIVI NNLO* and the NLO agree with the CDF data, 
as depicted by the gray (labeled NLO) and red (labeled 
NNLO*) band of Fig. 20. 

In this partial QCD NNLO calculation, some infrared 
divergences are not automatically regulated since the 
loop corrections at are missing. This imposed the 
introduction of an infrared cut-off, whose effect is be- 
lieved to decrease as pr increases, provided that the 
real-emissions dominates at large pr- This may not 
be exactly realised if double f-channel gluon exchanges 
have an important impact. In this case, the result could 
be overestimated. However, the agreement with the kr 
factorisation results gives us the hope that the infrared 
divergences are for a large part under control. 

This partial QCD NNLO calculation was also applied 
to the J /iff and the i//(2S) case [32] and the contribution 
was also found to be important. Yet, contrary to the Y 
case, there is still a small gap opening at large pr- At 
the leading order in v, the prediction for the J/tfr and the 
4^(28) are identical. The i/'(25) polarisation as function 
of pt is depicted by the red band (labeled NNLO*) band 
of Fig. 17. It is similar to that of NLO. 
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4.4. Polarization in photoproduction and at RHIC 
The J/t// polarization in photoproduction at HERA 

was also studied up to NLO in [124, 125]. The results 
showed that the transverse momentum, pr, and energy 
fraction, z, distributions of J/ 1// production do not agree 
well with the observations. The theoretical uncertain- 
ties on the z distributions of the J/iJ/ polarization pa- 
rameters in the target frame for various choices of the 
renormalization and factorization scales are too large to 
provide a definite prediction relative to the experimental 
data [126]. It is quite easy to understand that the uncer- 
tainty in the QCD effects could be rather large since the 
Pt of the J/t^ is quite small at HERA, I < pr < 5 GeV. 
Thus the theoretical prediction is not expected to agree 
well with the measurements. 

The NLO CSM polarisation of the J/i// produced in 
proton-proton collisions at RHIC at = 200 GeV has 
been studied [40]. The results show that the polarisation 
pattern at NLO in the helicity frame is in good agree- 
ment with the PHENIX data both in the centi-al [15] and 
the forward [12] regions when extra contribution from 
eg fusion and a data-driven range for the Xc contribu- 
tion to the y/i/r polarisation are considered. For the time 
being, nothing is known about f-wave feed downs at 
RHIC especially at large pr where xj starts to be large. 
This prevents any strong conclusion based on the com- 
parison with the data regarding the status of the CO con- 
tributions and a possible overestimation by the partial 
NNLO contribution. 

4.5. Perspectives on polarization studies 

To clarify the J/ifr polarisation puzzle, other observ- 
ables such as the measurements of J/tfr (T) hadropro- 
duction in association with a photon at the LHC can 
be considered. The NLO CS contribution was stud- 
ied [127] and the results show that the J/i// polariza- 
tion changes from transverse at leading-order to longi- 
tudinal at NLO. It is known that the CO contribution to 
the inclusive J /if/ hadroproduction may be one order of 
magnitude larger than the CS contribution, and the po- 
larization distribution may be dominated by the CO at 
NLO. In contrast, the CS contribution for J/t// + j pro- 
duction is at least of the same order as the CO contribu- 
tion. The polarization should then arise from both the 
CS and CO. Therefore, measurements of J/t/r produc- 
tion associated with a direct photon at hadron colliders 
could be an important test of the theoretical treatment 
of heavy-quarkonium production. The partial NNLO 
result [128], NNLO*, is about one order of magnitude 
larger than the NLO J/i// pr distribution. The predicted 
polarization are almost the same. Therefore, it will be 
reality check on the partial NNLO treatment. 



Moreover, to compare the polarisation measurements 
in different frames, we need theoretical predictions 
based on different polarisation frames such as the helic- 
ity frame (HX) and Collins-Soper frame. Along these 
lines, we would like to mention a recent pioner work on 
the polarization studies of T -H y in the kr factorization 
approach using different polarization frames [129]. 

Theoretical predictions of the new, proposed, frame- 
independent polarisation parameter A are needed. Of 
course, all the prediction should at NLO accuracy. Fur- 
thermore, the comparison between experimental data 
and theory must consider the experimental acceptance 
and efficiency. Experiments measure the net polarisa- 
tion of the specific cocktail of quarkonium events ac- 
cepted by the detector, trigger and analysis cuts. It will 
be very useful to have an event generator where the 
quarkonium decay into lepton pairs for inclusive pro- 
duction at NLO which could be embedded into Monte 
Carlo simulations. 

5. Open heavy-flavor production in pp and pA colli- 
sions 

Open heavy-flavor production in pp and pA colli- 
sions is tightly connected to the process of quarkonium 
(hidden flavor) production. Nevertheless, it is in it- 
self one of the most interesting processes to study from 
theoretical and experimental points of view (see, e.g., 
Refs. [130, 131] for a review on recent results). In 
fact, the production of heavy quarks allows us to test 
fundamental concepts, such as perturbative QCD, fac- 
torization and non perturbative power corrections. It 
also represents a mandatory benchmark for quarkonium 
production and parameters derived in open heavy-flavor 
production serve as the basis for calculations of quarko- 
nium production. 

Thanks to the massiveness of heavy quarks, the pro- 
duction cross section can be calculated analytically 
in perturbative QCD down to pr = at the par- 
tonic level [132]. However, differential distributions 
and event shapes exhibit large logarithmic contribu- 
tions, typically corresponding to soft or collinear par- 
ton radiation, which need to be resummed to all or- 
ders [133, 134, 135]. Furthermore, dead-cone effects 
suppress gluon radiation around the heavy-quark direc- 
tion, which has relevant phenomenological implications 
on both parton- and hadron-level spectra. As far as 
heavy-hadron production is concerned, one can model 
non perturbative effects by means of a fragmentation 
function, depending on a few tunable parameters, or 
by including power corrections in a 'frozen' or 'effec- 
tive' strong coupling constant [136, 137, 138]. Reliable 
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hadronization models should also be capable of distin- 
guishing the spin of the produced hadrons, e.g. sepa- 
rating the D* from the D, and describing both baryons 
and mesons and the multiplicity ratios A^/D or hi, IB. 
Typically, hadronization models, such as the cluster 
[139] or string [140] models, are implemented within 
the framework of Monte Carlo generators, such as HER- 
WIG [141], PYTHIA [142] orMC@NLO [143] and are 
tuned to heavy-hadron data from e^e" experiments, e.g. 
LEP and SLD data (see, for example, [144] for a tuning 
to /7-quark fragmentation). Multi-purpose generators 
have been available for several years for lepton/hadron 
collisions in the vacuum and have been lately modified 
to include the eff'ects of dense matter [145, 146] (see 
also [147, 148, 149, 150, 151] for generators specific 
for heavy-ion collisions). 

A profound knowledge of open heavy-quark produc- 
tion is crucial for the understanding of quarkonium phe- 
nomenology. The production, e.g., of the Jj^ in Non 
relativistic QCD (NRQCD) is described as the convolu- 
tion of the perturbative cc cross section and a non pertur- 
bative operator, possibly depending on the color of the 
cc state, according to the color-singlet and color-octet 
mechanisms (see e.g. [152]). 

In order to promote the formalism used to describe 
open heavy flavors to pA and ultimately AA collisions, 
one has to introduce medium-modification effects, tak- 
ing into account the existing dififerences between light 
and heavy quarks (hadrons) [153]. In particular, one of 
the most striking observations of heavy-ion collisions is 
the jet-quenching phenomenon, namely the suppression 
of hadron multiplicity at large transverse momentum 
{pr) with respect io pp processes [154, 155, 156, 157]. 

In the case of heavy hadrons, the measured suppres- 
sion is expected to be lower than that of light-flavor 
mesons, due to the dead-cone effect discussed above. 
This effect becomes less important in the high-pj- limit 
where pj » niQ. So far, suppression of heavy-flavored 
mesons in pA or AA collisions has only been mea- 
sured through leptons, i.e., the non-photonic electrons 
and muons from semi leptonic decays of heavy-flavored 
mesons [158, 159, 160]. Despite being a rather com- 
plex study, non-photonic electrons and muons have var- 
ious experimental advantages, the most important one 
being the possibility to deploy the fast triggers, in such 
a way to allow the experiments to collect large data sam- 
ples. On the other hand, they originate from both charm 
and bottom mesons with only few indirect experimental 
handles available to disentangle the two contributions. 

Suprisingly, the energy loss observed through non- 
photonic electrons is substantial akeady at moderate pj 
of ~3 GeV/c. Figure 21 shows this suppression via the 
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Figure 21: Nuclear modification factor (Raa) for non-photonic elec- 
trons originating from semileptonic decays of heavy-flavored mesons 
in the 0-10% most central AuAu collisions at RHIC measured by 
PHENIX [158]. This corresponds to the ratio of the invariant yields 
in AuAu and pp collisions normalized by the number of binary nu- 
cleon collisions. This measurement indicates a strong suppression of 
heavy-flavored mesons in central heavy-ion collisions. 



nuclear modification factor Raa defined as the ratio of 
invariant yields in AuAu over that in pp collisions times 
the number of binary nucleon collisions in the AuAu 
system at a given centrality. In the absence of any nu- 
clear or dense matter effects Raa is unity. Several mech- 
anisms to describe heavy-flavor energy loss in the dense 
medium have been so far proposed. Most prominent 
are collisional and radiative energy loss, typically calcu- 
lated in the weak-coupling regime. However, the mag- 
nitude of the observed suppression at RHIC is hard to 
accommodate in these models. Recent attempts, assum- 
ing strong coupling between the heavy-flavor mesons 
and the dense medium, and based on the AdS/CFT cor- 
respondence, lead to higher suppression factors, but fur- 
ther studies and experimental data are needed before 
any conclusions can be drawn. 

It is still an open question to what extent the energy- 
loss mechanisms, observed in open heavy-flavor pro- 
duction, affect the quarkonium yields in AA collisions. 
After including quark-antiquark recombination effects, 
the quenching of open heavy-flavored hadrons will also 
serve as a benchmark for quarkonia in heavy-ion colli- 
sions. 

Another crucial issue in heavy-flavor phenomenol- 
ogy concerns heavy-quark parton distributions in pro- 
tons and nuclei. In fact, processes with the produc- 
tion of heavy quarks are fundamental to constrain the 
PDFs. Higher-order corrections to the production of 
heavy quarks, with the implementation of soft/collinear 
resummation to improve the differential distributions. 
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will also be helpful to acquire information on heavy 
quark densities. The possible inclusion of medium ef- 
fects in the hard-scattering cross section, along with nu- 
clear parton distribution functions, will eventually lead 
to a prediction taking fully into account the modifica- 
tions induced by the dense matter. However, unlike col- 
lisions in the vacuum, there is no guarantee that factor- 
ization should also work in a medium. 

In the following, we shall review the main results dis- 
cussed within the Open Heavy-Flavor Working Group 
in the Quarkonium 2010 workshop. We shall discuss 
higher-order calculations of heavy-flavor production 
spectra at RHIC, taking particular care about the role 
played by the resummation on transverse-momentum 
distributions. We shall then present a few selected top- 
ics on specific features of heavy quarks and dead-cone 
effects, as well as novel developments on the modelling 
of energy loss from heavy quarks in dense matter and 
its possible relation with the jet-quenching observation. 
Later on, recent results on the production of photon 
plus heavy quarks, along with its relevance to constrain 
charm/bottom PDFs in proton and heavy ions, will be 
investigated. 

5.1. Open heavy-flavor production at RHIC and 
FONLL calculations 

As our understanding of heavy-ion collisions be- 
comes more mature, it becomes increasingly clear that 
heavy-flavor studies are a mandatory tool to further 
our understanding of the quark-gluon-plasma formed in 
high energy nuclear collisions at RHIC and now at LHC. 
Having reliable measurements of all aspects of heavy- 
flavor production is compelling. 

Experimentally, the total charm cross section is not 
well constrained. Measurements of charm meson pro- 
duction at the Tevatron could be conducted only at large 
transverse momenta {pj > 5 GeV/c) while measure- 
ments at lower energies exhibit rather large uncertain- 
ties. At RHIC, measurements of D mesons in con- 
junction with low-pj- non-photonic electrons and muons 
in dAu collisions [161] (STAR) indicate a rather large 
total charm cross section of ~1.3 mb. On the other 
hand, the non-photonic lepton measurements in pp col- 
lisions [158] (PHENIX) suggest a considerably smaller 
cross section of ~ 0.6 mb. In all cases assumptions 
about rapidity distributions and contributions from un- 
measured charm mesons and baryons (e.g. A^) have 
to be made, adding further contributions to the already 
large statistical and systematic uncertainties. Upgrades 
to the existing RHIC detectors (for an overview see 
[162]) will help to overcome these shortcomings. A 
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Figure 22: Comparison of total charm cross section measurements. 
The STAR and PHENIX results are given as cross section per binary 
collisions. Vertical lines reflect the statistical errors, horizontal bars 
indicate the systematic uncertainties (where available). From [162]. 

summary of current total charm cross sections (cr^j) is 
presented in Fig. 22. 

On the theory side, the state of the art for charm- 
quark production is NLO with the inclusion of «/ = 3 
active flavors [132]. A few ingredients of the NNLO 
corrections are available [163, 164], but the full NNLO 
calculation has not been completed yet. Another possi- 
ble approach to compute the cross section is to integrate 
the transverse momentum distribution computed in the 
so-called FONLL approximation, i.e. soft/coUinear 
resummation to next-to-leading logarithmic accuracy 
(NLL), matched to the NLO result. The total cross sec- 
tion is still NLO, but the charm quark is treated as an 
active flavor, and therefore nj - 4 in the calculation, 
which leads to a diff'erence in the result [135]. 

The FONLL calculation can of course be employed 
to predict differential distributions, in particular the 
transverse momentum one, wherein contributions ~ 
4 Xn^-^ipTlniQ) (LL) and ~ or* ln*^3(/?7-/'«e) (NLL), 
large for pj » mq, are summed up to all orders. The 
spectra of c-flavored hadrons are then obtained by con- 
voluting the charm distribution with the appropriate non 
perturbative fragmentation function, fitted, e.g., to LEP 
data. When doing the fit, the calculation of heavy-quark 
production in e^e^ annihilation must be carried out in 
the same perturbative approximation, i.e. NLO-hNLL, 
and the scales are to be consistently set [165, 166]. 

The main sources of uncertainty on such predictions 
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are the renormalization and factorization scales, the 
quark masses, entering in the perturbative calculation, 
parton distribution functions, such as the gluon density 
about the charm/bottom mass scale. As for the quark 
mass, it is also unclear whether in the calculation one 
should use the pole, the MS or even the hadron mass. 
For the purpose of the PDFs, the gluon distribution ex- 
hibits large uncertainties at small x, especially a strong 
dependence on the factorization scale. 
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Figure 23: Non-photonic electron measurements at RHIC (PHENIX) 
in pp collisions at = 200 GeV. (a) Invariant differential cross 
section of single electrons as a function of pj. (b) The ratio of 
FONLL/Data as a function of pj. The upper (lower) curve shows 
the theoretical upper (lower) limit of the FONLL calculation. From 
[158]. 

From the D- and B-hadron distribution one can de- 
rive the non-photonic electron spectrum at FONLL and 
compare it with RHIC data. Both the new PHENIX 
(shown in Fig. 23) and STAR pp data are compatible 
with the FONLL computation, within the theoretical 
and experimental uncertainties. However, the data fa- 
vors large cross sections, about 1.5 times the nominal 
FONLL calculation, just on the upper edge of the un- 
certainty band (Fig. 23b). Non-photonic electron spec- 
tra at low pt are notoriously difficult to analyze due to 
the increased photonic backgrounds, which increase the 
uncertainties for pj < 2 GeV/c considerably and make 
the evaluation of cr^g from electrons alone problematic. 

Furthermore, much work has been carried out to- 
wards the separation of charm from bottom contribu- 
tions to the lepton spectra in pp and pA collisions. Ex- 
perimental progress has been made using e-h and e-D 
correlations [167, 168] shown in Fig. 24. Again, results 
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Figure 24: Transverse momentum dependence of the relative contri- 
bution from B mesons to the non-photonic electron yields as mea- 
sured by PHENIX [167] and STAR [168]. Error bars are statistical 
and brackets/boxes are systematic uncertainties. The solid line is a 
FONLL prediction and the dotted lines represent the uncertainty on 
this FONLL prediction. 



are consistent with the FONLL calculation, although ex- 
perimental errors are currently still too large to allow an 
accurate separation of bottom and charm production in 
pp processes. 

The final answer has to come from direct measure- 
ments of charm through D mesons (if able to subtract 
the B ^ D contribution looking e.g. at the impact pa- 
rameter distribution) and that of bottom through B 
J/t// + X and/or /j-tagged jets or leptonic decays. These 
measurements will require high precision vertexing that 
will become available for both RHIC experiments in the 
near future [162]. 

5.2. Properties of gluon radiation off heavy quarks 

A heavy quark Q is characterized by the fact that 
its mass is much larger than the typical hadronization 
scale, say niQ » Aqcd, e-g-, in the MS renormalization 
scheme. The fact that the quark mass acts as a regulator 
of the coUinear singularity and that the strong coupling 
constant is evaluated at relatively large scales, depend- 
ing on the observable one is looking at, makes perturba- 
tive QCD applicable. 

When dealing with multiple emissions from heavy 
quarks, an essential difference with respect to light 
quarks is that gluon radiation is suppressed inside the 
forward angular cone with an opening angle mg/E, the 
so-called dead cone, with E being the heavy-quark en- 
ergy [169]. The dead cone implies the decrease of the 
overall energy loss and thus to the leading-particle ef- 
fect, i.e. heavy-hadron energy distributions are peaked 
at large x, with x being the energy fraction in the centre- 
of-mass frame [136]. In other words, after multi-parton 
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radiation and hadronization, a heavy hadron still car- 
ries a significant fraction of the parent-quark energy. 
Also, for X ^ I - Aqcd /inQ, one starts to be sensitive to 
non perturbative effects which can be included, e.g., by 
means of a fragmentation function or a frozen/effective 
coupling constant. Another relevant result, which can 
be related to the dead cone is that, unlike the naive 
parton-model expectation, the multiplicity difference 
between heavy and light hadrons is roughly independent 
of the centre-of-mass energy [170]. 

More generally, soft radiation obeys the Low- 
Barnett-Kroll (LBK) theorem [171, 172], according to 
which soft gluons are classical and therefore incapable 
of modifying the quantum numbers, such as the color, 
of a given system. Also, due to its classical nature, soft- 
parton emission is independent of the underlying hard- 
scattering process and of the quantum properties of the 
emitter. This implies that, if a process is forbidden by a 
given symmetry, such a veto cannot be relaxed by emit- 
ting extra soft photons/gluons. 

This property of soft radiation can be connected to 
the longstanding puzzle of J/t// production, whose mea- 
sured rate at large pr at the Tevatron was about 50 times 
larger than the theoretical prediction. In fact, the color- 
octet mechanism was formulated in order to explain 
such a discrepancy. In other words, such a model pre- 
dicts, at one loop in gluon-gluon fusion, the production 
of a J/tfr and a gluon, with the J/ 1// being in a color-octet 
state. 

However, as discussed above, in order to change the 
J/ iff color, the emitted gluon cannot be soft: according 
to the LBK theorem, the price to pay for the emission of 
a hard gluon is ~ (Aqcd/»Ic)^- 

Therefore, it is probably worthwhile investigat- 
ing alternative mechanisms which, besides the color- 
evaporation model and the color-octet model, can in- 
crease the yield of large-pj- J/t// production in hadron 
collisions and hence recover the agreement with the 
Tevatron data. In fact, the production of quarkonia at 
large transverse momentum has a very small cross sec- 
tion and therefore it is a rare configuration, wherein fluc- 
tuations are expected to play a role [173]. 

5.3. Heavy quarks and energy loss in dense matter 

A crucial issue in heavy-quark phenomenology is the 
description of energy loss in a dense medium and its re- 
lation to the quarkonium spectrum in matter. It is com- 
monly assumed that partons at large transverse momen- 
tum traversing the medium lose energy through incoher- 
ent radiation of gluons as well as collisional energy loss. 
The approach presented in [174] assumes that a fast par- 
ton in a medium undergoes multiple collisions, charac- 



terized by a moderate momentum transfer and large val- 
ues of the strong coupling constant. As for the coupling, 
non perturbative effects are taken into account by using 
a frozen coupling constant [136, 137]. Indeed, such a 
model manages to reproduce quite well the Raa ratio 
measured by the PHENIX collaboration, for AuAu col- 
lisions at 200 GeV over the full transverse-momentum 
range, up to a normalization /T-factor of about 2. 

As for the radiative energy loss, for light quarks 
one usually relies on the Baier-Dokshitzer-Mueller- 
Peigne-Schiff (BDMPS) approximation [175] which 
assumes static scattering centres, independent soft emis- 
sions and hadronization outside the medium. For heavy 
quarks one can rely to the Gunion-Bertsch approach 
[176]: the medium-induced gluon-radiation spectrum is 
corrected by means of terms depending on the heavy- 
quark mass and on the fictitious gluon thermal mass. 
Moreover, gluons radiated by massive quarks are re- 
solved in a smaller time with respect to light quarks, 
which implies a milder dependence on coherence effects 
and evolution variables for multiple radiation. There- 
fore, as a first approximation, one can even think of 
neglecting coherence effects when dealing with emis- 
sions from heavy quarks. The result is that if one adds 
both collisional and radiative energy loss contributions, 
a scaling factor of 0.6 is enough to reproduce the R^a 
data at RHIC. Future data from the LHC will possibly 
shed light on the issue of the energy loss; to this goal, it 
will be very useful having more exclusive probes, such 
as azimuthal correlations or tagged c- or ^-flavored jets. 
Another open issue concerns the gluon thermal mass 
and whether it can be theoretically predicted or it should 
rather be fitted to experimental data, like a non pertur- 
bative parameter 

5.4. Heavy quark plus photon production as a probe of 
the heavy quark PDF 

The production of a direct photon accompanying a 
heavy quark at hadron colliders is an important process, 
also relevant to constrain the heavy-quark density. As 
it escapes confinement, a photon can be exploited to 
tag the hard scattering; also, its transverse momentum 
distribution gives meaningful information on the heavy- 
quark production process and PDF. 

At leading order, i.e. 0{aas), there is only one hard 
scattering process, namely gQ — > Qy. At NLO in the 
strong coupling constant, i.e. Oiaa^^), several subpro- 
cesses contribute to Qy production, with gluons, heavy 
and light quarks in the initial state [177]. Considering, 
e.g., charm production, such processes will help to shed 
light on the charm distribution and whether it is radia- 
tively generated by means of gluon emissions, in which 
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case c(x, Q^) ^ g{x, Q^), or there is an intrinsic charm 
contribution to the proton. In fact, there exist light cone 
models, wherein intrinsic charm is relevant at large x 
[178], and sea-like models, with c(x) ^ u(x) + d(x). 
As for the comparison between the photon spectra pro- 
duced in association with charm and bottom quarks, 
there are remarkable differences both at LO and at NLO. 
At large transverse momenta, the discrepancy between 
LO and NLO gets larger [177]. 

The NLO calculation has been compared with DO 
data [179] on the photon transverse momentum in Qj 
events: the agreement is pretty good in the case of 
by production, whereas discrepancies are present for 
charm production, with the NLO computation yielding 
a lower cross section for moderate and large values of 
Pt- Such disagreement can be traced back to the possi- 
ble existence of intrinsic charm in the proton. Using the 
Brodsky-Hoyer-Peterson-Sakai (BHPS) PDF set [178] 
slightly improves the comparison, although meaningful 
differences are still present at very large px- The dis- 
crepancy is instead milder when comparing the ratio 
of charm and bottom cross sections. At the LHC, one 
can still investigate the presence of the intrinsic charm. 
However, since the probed values of x are on average 
smaller than at the Tevatron, one is basically sensitive 
to the sea-like model. Also, the forward-rapidity region 
will be particularly suitable to discriminate between ra- 
diative and intrinsic contributions to the charm PDF. 

The LO/NLO calculations for b/c + y production can 
be extended to pA collisions, e.g. proton-lead processes 
at the LHC. A first approximation consists in using the 
same amplitudes as in the vacuum and convoluting them 
with nuclear parton densities, such as the EPS [55], 
HKN [180] and nCTEQ [181] sets. UnUke the pro- 
ton case, the gluon distribution in a nucleus is poorly 
known, and it is one of the main differences among the 
nuclear PDF sets. In fact, possible measurements at the 
LHC of 7 + c/b final states will help to investigate both 
gluon and charm/bottom distributions in dense matter 

Likewise, similar studies can be carried out for 
deuterium-gold collisions at RHIC. Since the Bjorken- 
X regions probed at RHIC and LHC are complementary, 
these measurements will help to discriminate among the 
nuclear parton distribution sets. As for the hard scat- 
tering, how to implement medium modifications, for 
both light and heavy quark production, is still an open 
issue; moreover, collinear factorization is an assump- 
tion that has not yet been proved. The BDMPS and 
Gunion-Bertsch approximations discussed above work 
well only for soft/collinear emissions. Indeed, there 
are prescriptions for hard/large-angle radiation, but a 
thorough implementation and phenomenological analy- 



sis on medium-modified hard scatterings is still missing. 
Convoluting possible medium-induced hard matrix ele- 
ments with nuclear parton densities will ultimately lead 
to cross sections fully including dense-matter effects. 

6. Quarkonia as a Tool 

6.1. Proton-Proton Collision Studies 

Since the start of LHC operation, new experimental 
data on quarkonium production has revived the interest 
in understanding the production mechanisms of various 
quarkonium states in hadronic collisions. All the mod- 
els describing various sources of onia in hadronic colli- 
sions share the common basis: the factorisation between 
the hard collison subprocess and the parton-parton col- 
lision luminosity, calculated as a convolution of the par- 
ton distribution functions (PDFs). 

Clearly, the range of accessible xi^2 depends on the 
rapidity interval covered by the experiments. For char- 
monium at the Tevatron, with a partonic center-of-mass 
energy VI - 3.5 GeV, sfs - 2 TeV, the CDF range is : 

yr^O ^ Jt:i,2 ^ 1.8 X 10"^ 

while for DO: 

-1-6<3'<1.6 => Jt;i,2 ^ (0.36 - 8.9) X 10"^ 

Their values of xi 2 lie within the range where the PDFs 
are well-established. Hence the uncertainties on the par- 
ton luminosities at the Tevatron are fairly small. 

At the LHC, not only is the energy higher by a factor 
of 3.5, but the rapidity coverage of the various experi- 
ments is also significantly wider, especially if one com- 
bines the results from ATLAS and CMS with those of 
LHCb and ALICE. One get for ATLAS and CMS: 

- 2.4 < y < 2.4 ^ xi=^ (0.04 - 6.0) x 10"^ 
X2 ^ (6.0 - 0.04) X 10"^ 

while for LHCb: 

2 < y < 5 ^ xi ^ (4 - 80) ■ 10"^ 

X2 ^ (0.07 - 0.003) X 10"^ 

ALICE has complementary coverages {\y\ < 1 & 2.5. < 
y < 4.0). For the kinematics of LHC experiments, 
the accessible x range of charmonium production mea- 
surements is up to three orders of magnitude below the 
"safe" value of lO^^*. Processes usually used for de- 
termining gluon PDFs, such as jets, prompt photons or 
open charm pairs, are for various reasons unlikely to be 
useful at these extremely low x, leaving charmonium as 
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the only source of information on gluon PDFs in this 
area. 

The smallest value for bottomonium is about 
10 GeV, three times larger than for charmonium. This 
restricts the "reach" of T studies to xi_2 values three 
times larger than those contributing to Jj^. However, 
this is still significantly smaller than the currently acces- 
sible range, and thus should be useful as a consistency 
check and maybe also to test the evolution (provid- 
ing data with a higher for similar xb, which may be 
helpful to differentiate between a BFKL- and DGLAP- 
like evolution). 

The main experimental difficulties include a reduced 
acceptance for J 1 41 when both pj and |>i are small, 
and the necessity of suppressing non-prompt charmo- 
nia (which come from larger s). Theoretical complica- 
tions include proper and consistent treatment of intrinsic 
parton pj, initial-state radiation, and higher-order per- 
turbative effects, as well as the possibility of gluon re- 
combination at high densities. So concerted efforts from 
theorists and experimentalists will be needed to extend 
our measurements of the PDFs to much smaller x than 
prior to the LHC. This is necessary for understanding 
the dynamics of quarkonium hadronic production. 

6.2. Nuclear Collision Studies 

From statistical QCD we know that strongly inter- 
acting matter undergoes a deconfinement transition to 
a new state, the quark-gluon plasma (QGP). The objec- 
tive of high energy nuclear colhsions is to produce and 
study the QGP under controlled conditions in the lab- 
oratory. How then can we probe the QGP - what phe- 
nomena provide us with information about its thermo- 
dynamic state? Here the ultimate aim must be to carry 
out ab initio calculations of the in-medium behaviour 
of the probe in finite temperature QCD. Quarkonia may 
well provide the best tool for this presently known, but 
in nuclear collisions other phenomena wiU also lead to 
modified quarkonium production. As a consequence, 
parton distribution changes, parton energy loss and cold 
nuclear matter effects (absorption) on quarkoniirai pro- 
duction must be accounted for before any QGP studies 
become meaningful. For the sake of definiteness, in this 
section we shall refer to the non-QGP effects as normal, 
in contrast to the anomalous suppression (or enhance- 
ment) we want to look for. 

Any modifications observed when comparing 
J I if/ production in AA collisions to that in pp inter- 
actions thus have two distinct origins. Of primary 
interest is obviously the effect of the secondary medium 
produced in the collision - this is the candidate for 
the QGP we want to study. In addition, however, the 



presence of the "normal" initial and final state effects 
will presumably also affect the production process and 
the measured rates. This ambiguity in the origin of any 
observed Jlij/ suppression will thus have to be resolved. 

A second empirical feature to be noted is that the 
observed y/i/r production consists of directly produced 
\S states as well as of decay products fvomXcO^) and 
i/''(2S) production. The presence of a hot QGP affects 
the higher excited quarkonium states much sooner (at 
lower temperatures) than the ground states. This results 
in another ambiguity in observed Jjij/ production - are 
only the higher excited states affected, or do all states 
suffer? 

An ideal solution of these problems would be to mea- 
sure separately Jlil/,Xc and if/' production first inpA (or 
dA) collisions, to determine the effects of cold nuclear 
matter, and then measure, again separately, the produc- 
tion of the different states in AA collisions as function 
of centrality at different collision energies. 

Until such studies become available, we resort to a 
more operational approach, whose basic features are: 

• We assume that the Jlij/ feed-down rates in pA and 
AA are the same as in pp, i.e., 60 % direct Jjij/ , 30 
% decay of x(lf'), and 10 % decay of il/'{2S). 

• We specify the effects due to cold nuclear matter 
by a Glauber analysis of pA or dA experiments in 
terms of cr'^^^^ for i =JI^, Xc, This a'^^ is not 
meant to be a real cross section for charmonium 
absorption by nucleons in a nucleus; it is rather 
used to parametrize all initial and final state nuclear 
effects, including shadowing/antishadowing, par- 
ton energy loss and pre-resonance/resonance ab- 
sorption. 

• In the analysis of AA collisions, we then use o-'^^^^ 
in a Glauber analysis to obtain the predicted form 
of normal J/tfr suppression. This allows us to iden- 
tify anomalous J/iJ/ suppression as the difference 
between the observed production distribution and 
that expected from only normal suppression. We 
parametrize the anomalous suppression through 
the survival probability 

JANildy)^ ^^^^ 

idNi/dy)GUnbeT 

for each quarkonium state i. 

With the effects of cold nuclear matter thus accounted 
for, what form do we expect for anomalous J/i// sup- 
pression? If AA colhsions indeed produce a fully equi- 
librated QGP, we should observe a sequential suppres- 
sion pattern for J /if/ and Y, with thresholds predicted 
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(quantitatively in temperature or energy density) by fi- 
nite temperature QCD [183, 184, 185, 186]. The result- 
ing pattern for the is illustrated in Fig. 25. 
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Figure 25: Illustration of sequential Jlijj suppression. 

Its consequences are quite clear. If, as present sta- 
tistical QCD studies indicate, the direct J 1x1/ survives 
up to about 2 Tc and hence to e > 25 GeV/fm-', then all 
anomalous suppression observed at SPS and RHIC must 
be due to the dissociation of the higher excited states 
Xc and i/r' . The suppression onset for these is predicted 
to lie around e - I GeV/fm^, and once they are gone, 
only the unaffected J I if/ production remains. Hence the 
J I ill survival probability for central AuAu collisions at 
RHIC should be the same as for central PbPb collisions 
at the SPS. 

A further check to verify that the observed y/i/r pro- 
duction in central collisions is indeed due to the un- 
modified survival of the directly produced 15 state is 
provided by its transverse momentum behaviour Initial 
state parton scattering causes a broadening of the pj dis- 
tributions of charmonia [187, 188, 189]: the gluon from 
the proton projectile in pA collisions can scatter a num- 
ber of times in the target nucleus before fusing with a 
target gluon to produce a cc. Assuming the gluon from 
the proton to undergo a random walk through the target 
leads to 

{pl)pA^{pl)pp+N^5a (15) 

for the average squared transverse momentum of the ob- 
served y/i/r. Here specifies the number of collisions 
of the gluon before the parton fusion to cc, and 5o the 
kick it receives at each collision. The collision number 
A^^ can be calculated in the Glauber formalism; here 
(Tabs has to be included to take the presence of cold nu- 
clear matter into account, which through a reduction of 
J 1x1/ production shifts the effective fusion point further 
"down-stream" [190]. 

In AA collisions, initial state parton scattering occurs 
in both target and projectile and the corresponding ran- 



dom walk form becomes 

{p]-)AA^{p]-)pp+N^5Q, (16) 

where now N^^ denotes the sum of the number of col- 
lisions in the target and the projectile prior to parton fu- 
sion. It can again be calculated in the Glauber scheme 
including cr^<o^. The crucial point now is that if the ob- 
served J I ill's, in central AA collisions are due to undis- 
turbed IS production, then the centrality dependence of 
the pt broadening is fully predicted by such initial state 
parton scattering [186]. In contrast, any onset of anoma- 
lous suppression of the y/i/r would lead to a modifica- 
tion of the random walk form [190]. 

In Fig. 26 we summarize the predictions for Jlxjj sur- 
vival and transverse momentum behaviour in AA col- 
lisions at SPS and RHIC, as they emerge from our 
present state of knowledge of statistical QCD. Included 
are some preUminary and some final data; for a discus- 
sion of the data analysis and selection, see Ref. [186]. 
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Figure 26: Jj\p survival and transverse momentum at SPS and RHIC. 

We conclude that the present experimental results are 
compatible with the present information from statistical 
QCD. This was not the case previously, and such a con- 
clusion can be drawn today because of several changes 
in our theoretical and experimental understanding. 
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• Finite T lattice QCD suggests direct suppres- 
sion at energy densities beyond the RHIC range; 
previous onset values were much lower (see, e.g., 
Ref. [182]). 

• SPS Inin data suggest an onset of anomalous sup- 
pression at e ^ 1 GeV/fm-'; previous onset values 
from PbPb and SU interactions were considerably 
higher, with e ^ 2 - 2.5 GeV/fm^ (see, e.g., Ref. 
[191]). 

• Within statistics, there is no further drop of the 
Jlif/ survival rate below 50 - 60 %, neither at RHIC 
nor at the SPS; a second drop in very central SPS 
PbPb data (see, e.g., Ref. [191]) is no longer con- 
sidered viable. 

6.2. 1. J I tji Enhancement by Regeneration 

In this section we want to consider an alternative ap- 
proach which can so far also account for the available 
data. The basic idea here is that the medium produced 
in nuclear collisions is not a QGP in full equilibrium but 
rather one which is oversaturated in its charm content. 

A crucial aspect in the QGP argumentation of the 
previous sections was that charmonia, once dissociated, 
cannot be recreated at the hadronization stage since the 
abundance of charm quarks in an equilibrium QGP is 
far too low to allow this. The thermal charm quark pro- 
duction rate, relative to that of light quarks, is 

— ^ exp{-m,/r,) ^ 6 X 10"^ (17) 
qq 

with nic - 1.3 GeV for the charm quark mass and 
Tc - 0.175 GeV for the transition temperature. The 
initial charm production in high energy nuclear inter- 
actions, however, is a hard non-thermal process, and 
the resulting c/c production rates grow with the number 
(A^coii) of nucleon-nucleon collisions. In contrast, the 
light quark production rate grows (at least in the present 
energy regime) essentially as the number (A^pan) of par- 
ticipant nucleons, i.e., much slower The initial charm 
abundance in AA collisions is thus much higher than the 
thermal value; we illustrate this in Fig. 27 for A - 200 
as function of the collision energy ^[s. What happens 
to this excess in the course of the collision evolution? 

The basic assumption of the regeneration approach 
[192, 193, 194] is that the initial charm excess is main- 
tained throughout the subsequent evolution, i.e., that the 
initial chemical non equilibrium will persist up to the 
hadronization point. If that is the case, a c from a given 
nucleon-nucleon collision can, at hadronization, com- 
bine with a c from a different collision ("off-diagonal" 
pairs) to recreate a Jfij/. This pairing provides a new sec- 
ondary statistical charmonium production mechanism. 
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Figure 27: Thermal vs. hard charm production. 

in which the c and the c in a charmonium state have 
different parents, in contrast to the primary dynamical 
production in a pp collision. At sufficiently high energy, 
this mechanism will lead to enhanced Jjip production 
in AA collisions in comparison to the scaled pp rates. 
When should this enhancement set in? 

In recent work [192, 193, 194], it is generally as- 
sumed that the direct Jlx}/ production is strongly sup- 
pressed for e > 3 GeV/fm-'. This is evidently in contrast 
to the statistical QCD results discussed in the previous 
sections; however, we recall the caveat that the tempera- 
ture dependence of the charmonium widths is so far not 
known. Moreover, it is of course always possible that 
the medium produced in nuclear collisions is quite dif- 
ferent from the quark-gluon plasma of statistical QCD. 

Next, it is assumed that the regeneration rate is de- 
termined by statistical combination in a QGP in kinetic 
equilibrium, with or without wave function corrections. 
In other words, if a c and a c meet under the right kine- 
matic conditions, they are taken to form a Jjij/. An evo- 
lution towards a QGP in chemical equilibrium would 
also allow annihilation at this point. 

To account for the y/t/f production rates observed at 
RHIC, it is then assumed that the new statistical produc- 
tion just compensates the proposed decrease of the di- 
rect primary 15 production, as illustrated in Fig. 28. At 
the LHC, with much higher energy densities, one should 
then observe a y/i^ enhancement relative to the rates ex- 
pected from scaled pp results. 

We thus have to find ways of distinguishing between 
the two scenarios discussed here: the sequential sup- 
pression predicted by an equilibrium QGP or a stronger 
direct suppression followed hy a J/tfr regeneration in a 
medium with excess charm. Fortunately the basic pro- 
duction patterns in the two cases are very different, so 
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Figure 28: Illustration of regeneration of Jjtfj production. 

that one may hope for an eventual resolution. 

The overall survival probability in the two cases 
is illustrated in Fig. 29. Sequential suppression provides 
a step-wise reduction: first the higher excited charmo- 
nium states are dissociated and thus their feed-down 
contribution disappears; at much higher temperature, 
the y/i/r itself is suppressed. Both onsets are in principle 
predicted by lattice QCD calculations. In the regenera- 
tion scenario, the thermal dissociation of all "diagonal" 
7/ (/r production is obtained by extrapolating SPS data to 
higher energy densities. The main prediction of the ap- 
proach is therefore the increase of 7/ ijj production with 
increasing energy density. Ideally, the predictions for 
the LHC are opposite extremes, providing of course that 
here the feed-down from B-decay is properly accounted 
for. 
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Figure 29: Jjijj survival: sequential suppression vs. regeneration. 

The expected transverse momentum behaviour in the 
two cases is also quite different. In the region of full 
y/i/f survival, sequential suppression predicts the nor- 
mal random walk pattern specified through pA stud- 
ies; the eventual dissociation of direct J I if/ states then 
leads to an anomalous suppression also in the aver- 
age [190]. Regeneration alone basically removes 
the centrality dependence since the different partners 



come from different collisions. It is possible to intro- 
duce some small centrality dependence [195], but the 
random walk increase is essentially removed. The re- 
sulting behaviour is schematically illustrated in Fig. 30. 
More generally, the quarkonium momentum distribu- 
tions, whether transverse or longitudinal, should in the 
regeneration scenario be simply a convolution of the 
corresponding open charm distributions; this provides 
a further check [195]. 
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Figure 30: /)7--behaviour: sequential suppression vs. regeneration. 

6.2.2. Conclusions 

• In statistical QCD, the spectral analysis of quarko- 
nia provides a well-defined way to determine tem- 
perature and energy density of the QGP. 

• If nuclear collisions produce a quark-gluon plasma 
in equilibrium, the study of quarkonium produc- 
tion can provide a direct way to connect experi- 
ment and statistical QCD. 

• For a QGP with surviving charm excess, off- 
diagonal quarkonium formation by statistical com- 
bination may destroy this connection and instead 
result in enhanced y/f/r production. 

7. New observables in quarkonium production 

In order to better understand the mechanisms respon- 
sible for heavy quarkonium production, it is now essen- 
tial to introduce new observables. The theoretical un- 
certainties on the predicted yields are too large to draw 
any strong conclusions on production models. More- 
over, polarization studies are known to be experimen- 
tally challenging and not always possible over the full 
phase-space region even if the yield is otherwise mea- 
sured. The LHC will certainly provide new information 
on the yields and polarizations. In addition, the LHC 
is well positioned to make unique new measurements 
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previously unavailable due to unsufRcient luminosity, 
center-of-mass energy, detector resolution or even man- 
power. 

7.1. Hadronic activity around quarkonium 

The first observable that we discuss, hadronic activ- 
ity around quarkonium [196], is, in fact, not completely 
new. Indeed, UAl compared their charged-track dis- 
tributions with Monte Carlo simulations of a J/t// com- 
ing from a B decay and a J/iff coming from a Xc de- 
cay [197, 198]. In the early nineties -before the Teva- 
tron results- Xc feed-down was still expected to be the 
major source of prompt J/tfr production. However, if 
we assume either that production proceeds through CO 
transitions or through CS transitions at higher-orders, 
we now expect more complex distributions, even for the 
prompt yield. It is therefore not clear how accurately 
such methods - recently revived by STAR in [13] - can 
determine the B-feeddown to // tfr, other than with mea- 
surements of displaced vertices typical of B decays. In- 
stead, the hadronic activity around quarkonium could 
be a good discriminant between CO and CS contribu- 
tions to quarkonium production since gluon emission 
during CO transitions would increase the hadronic ac- 
tivity. This study is, however, difficult to carry out in 
practice [196]. 

7.2. Associated production 

We therefore urgently need more observables easier 
to predict that test the many production models avail- 
able [107, 109]. We argue that the study of associated 
production channels, ijj + cc and T + bb, first in pp 
collisions, then in pA and AA, fills both these require- 
ments. These reactions are potentially very interesting 
since the LO description already shows leading pr be- 
haviour. We expect that higher-QCD corrections will 
be less important and the cross section will thus be less 
sensitive to the renormalisation scale, allowing for more 
precise predictions. 

Further motivation for such studies is the amazingly 
large fraction of J/tfr produced in association with an- 
other cc pair at B-factories. For instance, the Belle col- 
laboration first found [199] 



cr{e^e —» 7/^^-1- cc) 
o-ie+e- Jlip + X) 



= 0.59;°:[^±0.12. (18) 



More recently, the Belle collaboration reported new 
measurements [37]: 

o-(J/4r + X) = 1.17 + 0.02 + 0.07 pb, 
o- (J/i// + cc) = 0.74 + 0.08!^:!^^ pb , (19) 
cr ( J /if, + Xnon cc) = 0.43 ± 0.09 ± 0.09) pb. 



which confirm that associated J /if/ + cc production is the 
dominant production mechanism at B factories. 

Until now, it is unknown whether such a high associ- 
ated fraction also holds for hadroproduction. Analyses 
at the Tevatron (CDF and DO) and at RHIC (PHENIX 
and STAR) are feasible, in addition to the LHC. The 
predicted Tevatron Run 11 integrated cross sections at 
= 1-96 TeV are significant [23] : 



cr(7/(/r + cc) X Br(£^r) ^ 1 nb, 
cr(T + bb) X Br(tr) ^ 1 pb. 



(20) 



Neglecting the likely reduction of the CO matrix ele- 
ments needed to comply with the e^e^ results and in- 
duced by higher-order QCD corrections to the color- 
singlet contributions (see section 2), the integrated 
cross sections for associated J/tfr + cc production were 
found [200] to be dominated by the CS part. This is also 
the case for the pr distribution up to at least 5 GeV for 
tfr and 10 GeV for T. This clearly means that such ob- 
servables directly probe the CS mechanism. For most of 
the other observables, it is not as clear whether the CO 
contribution dominante over the CS one and they cannot 
be regarderd as a clean probe of CS constributions due 
to the inherent uncertainty from the CO yield. 

If the effect of CO transitions is indeed negligible 
for inclusive T production, T's produced in association 
with a bb pair would be unpolarised at the LHC, inde- 
pendent of pt- 

Besides discriminating between the CO and the CS 
transitions, the i/r yield in association with a cc should 
show an a priori completely different sensitivity to the 
Xc feed-down than the inclusive yield. The same holds 
true for Y + bb relative toxb feed-down. The CS f -state 
yield is expected to be smaller than the S -state one since 
they are suppressed by powers of the relative velocity v. 
Contrary to the inclusive case, no additional gluon is 
needed to attach to the heavy quark loop to produce the 
i// (or Y) compared to P-states. 

For the CO transitions, associated production of Xc + 
cc can occur via gg — > gg where both final-state gluons 
split into cc pairs. One of them then hadronizes into 
axc via a CO transitions. This contribution is certainly 
suppressed up to pr - 20 GeV. At larger pr, a dedicated 
calculation is needed. However, this mechanism could 
be easily disentangled from the CS contributions since 
both c quarks are necessarily emitted back to back from 
the Xc and thus to the J/ tf/. 

The non-prompt signal would originate as usual from 
gg — > bb, where one b quark decays into a ifr + X. Usu- 
ally, hadronisation of the b produces the i/r with light 
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quarks only. This means that we have a single c quark 
in the event, produced from the decay of the recoiling b 
quark and therefore back to back from the i//. The non- 
prompt signal could then be simply suppressed down 
by searching for a D meson near the i/'. It is possible 
that hadronisation of the b produces both the ip and a 
D meson. In this case, kinematic cuts would not help 
suppress the non-prompt yield. Fortunately, this final 
state is a priori suppressed compared to the first case 
and even further suppressed relative to the direct yield. 
Indeed, there is no gain in the pr dependence since both 
the gg — » bb and gg ^ i// + cc cross sections scale like 
Pj.*. A cross check sizing up the non-prompt yield with 
a displaced vertex measurement would be instructive. 

Associated production has also been studied in direct 
77 collisions in ultraperipheral collisions (UPC) [201]. 
At least for direct 77 collisions, associated production is 
the dominant contribution to the inclusive rate for pr > 
2 GeV/c. 

To conclude, studies can be carried on by detecting 
either the "near" or "away" heavy-quark with respect 
to the quarkonia. There are, of course, different ways 
to detect the D, B, or b-jet, ranging from a displaced 
vertices to detection of semileptonic decays to e or p. 

As mentionned in section 4, associated production of 
J/i// or Y with a prompt photon may also be investigated 
at the LHC. Interestingly, contrary to the inclusive case, 
the C - +1 CO should impact the large pj region, while 
C = -1 CO the small and mid pr region. The NLO CS 
contribution was recently studied [127] and, contrary to 
the inclusive case, it was shown that for this process the 
CO and CS induced yield should be similar, even with 
the largest possible CO matrix elements. At NNLO, the 
CS yield is expected to be further enhanced [128] and 
to dominate over the CO yield. The measurements of 
J/i// and T production associated with a direct photon at 
hadron colliders could thus be an important test of the 
theoretical treatment of heavy-quarkonium production. 

The background to these processes should also be 
taken into account. The Quarkonium event genera- 
tor Madonia [202], imbedded in Madgraph [203], will 
surely be of a great help in achieving this. In any case, 
we hope that such measurements would provide clear 
information on the mechanisms at work in quarkonium 
production. 

8. Conclusion 

Quarkonium production measurements in pp colli- 
sions have motivated many theoretical investigations. 
However, despite advances, there is still no clear picture 
of the quarkonium hadroproduction mechanism. Such 



a mechanism should be able to explain both the cross 
section and the polarization measurements at the Teva- 
tron and RHIC, as well as comparable measurements at 
B factories and lepton-hadron colliders. Presently, we 
have several different approaches at our disposal, such 
as the Color Singlet Model (CSM), in which the QQ 
is produced color neutral at short distances, and Non- 
Relativistic QCD (NRQCD), where quarkonium pro- 
duction can also proceed via creation of color-octet QQ 
pairs, via the Color Octet Mechanism (COM). These ap- 
proaches have both advantages and disadvantages. 

NRQCD has successfully explained the excess in 
charmonium production reported by the CDF Collab- 
oration, orders of magnitude larger than the LO CSM. 
However so far it fails to provide a consistent descrip- 
tion of production and polarization in pp collisions and 
the rates at B factories. 

Higher order corrections in ffj have only recently 
been calculated. The CS corrections at high pr are very 
large given that the corrections open new production 
channels which are relevant at high pr- On the contrary, 
CO calculations show that the cross sections do not in- 
crease much when these corrections are taken into ac- 
count since NLO CO contributions do not open any new 
channels. The contributions of the NLO CS corrections 
reduce the discrepancy between the CSM cross sections 
and the CDF data, but still fall too steeply at high pr to 
describe this region successfully. The NNLO* contribu- 
tion may be able to fill the gap between calculations and 
high pt data. At RHIC, the high pr STAR data seems 
to favor the COM over the CSM, while at low pr the 
PHENIX data are well reproduced by both approaches. 
It is, however, worth noting that little is known about 
feed down contributions at RHIC, while there does not 
exist any NLO theoretical predictions for the Xc yield, 
in either approach. 

In addition, the LO NRQCD calculation predicts a 
sizable transverse polarization rate for large pr J I if/ 
whereas the Tevatron CDF measurement at Fermilab 
displays a slight longitudinal polarization at large pj 
for the prompt yield. This problem is still present 
in calculations at NLO, since these corrections only 
slightly change the pj dependence of Jlijj produc- 
tion rate and polarization with respect to the LO 
calculations. However, the CSM results drastically 
change from transverse-polarization dominance at LO 
to longitudinal-polarization dominance at NLO, nearly 
in agreement with the data on 7/ i/r polarization produced 
in pp collisions both at the Tevatron and RHIC, if the 
polarization of the fmmxc feed down is assumed. 

In order to solve this puzzle, improved measurements 
are needed. So far, most experiments have presented 
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results based on a fraction of the physical information 

derivable from the data: only one polarization frame 
is used and only the polar projection of the decay an- 
gular distribution is studied. These incomplete results 
prevent model-independent physical conclusions. In the 
forthcoming LHC measurements, it is important to ap- 
proach the polarization measurement as a multidimen- 
sional problem, determining the full angular distribution 
in more than one frame. 

The feed downs from more massive states, such as 
higher-mass quarkonium states, cannot yet be taken into 
account in most of the calculations. While there are 
CDF measurements of direct production, the Xc feed 
down can significantly influence the polarization of the 
prompt J/ifr yield. A better way to avoid the feed-down 
problem is to study i//' and T(3S) where no important 
feed down from higher excited states exist. 

The behavior of quarkonium production in pA col- 
lisions due to Cold Nuclear Matter (CNM) effects is 
not yet understood. Puzzling features in proton-nucleus 
data put forward new aspects of charmonium physics in 
nuclear reactions, namely the role of cold nuclear matter 
effects (CNM). Two CNM effects are considered to be 
of particular importance: the modification of the initial 
parton disttibutions (PDFs) due to the nuclear environ- 
ment, an initial- state effect known as shadowing; and 
the breakup of cc pairs after multiple scatterings with 
the remnants of target nucleus, referred to as the nuclear 
absorption. 

Parton shadowing in the target nucleus may suppress 
(or enhance, in case of antishadowing) the probability 
of producing a J/i//. Its effect depends on the kinemat- 
ics of the cc pair production. The theoretical expecta- 
tions for gluon shadowing are quite diverse. A weak 
shadowing effect was predicted by the dipole approach 
and in some analysis of DIS data such as the one de- 
veloped by deFlorian and Sassot at leading and next- 
to-leading order - the nDS and nDSg parameterizations 
-. However, strong gluon shadowing and antishadowing 
was obtained in the EKS98, EPS08, and EPS09 param- 
eterizations by Eskola and collaborators. We note that 
some pA data were included in the EPS08 and EPS09 
analyses. The nDSg and EKS98 parameterizations are 
compatible for x < 10"^. 

Secondly, there is the effect of the nuclear absorp- 
tion. The strength of the interaction of the evolving cc 
pair with the target nucleons can break up the pair and 
consequently suppress the J/t// yield. The intensity of 
this effect may also depend on the quantum states of the 
cc pair at the production level (color octet or color sin- 
glet), and on the kinematic variables of the pair. Indeed, 
a colorless cc pair created in a hard reaction is not yet 



charmonium since it does not have a fixed mass. It takes 

time to disentangle ground state charmonium from its 
excitations, tf ~ l/(w^, - m^/^)- This time scale is in- 
terpreted as the formation time of the charmonium wave 
function. At low energies, this time is shorter than the 
mean nucleon spacing in a nucleus so that the forma- 
tion process can be treated as instantaneous. At much 
higher energies, the formation time is long and the ini- 
tial size of the cc dipole is "frozen" by Lorentz time 
dilation while propagating through the nucleus. 

According to the above arguments, at high energy, 
the heavy state in the projectile should undergo coherent 
scattering off the nucleons of the target nucleus, in con- 
trast to the incoherent, longitudinally ordered scattering 
that takes place at low energies. This should lead to a 
decrease of cr^^^ with increasing energy -^saw- Compi- 
lation and systematic study of many experimental data 
indicate that crJJ^^ appears either constant or decreasing 
with energy, following the most recent theoretical ex- 
pectations. 

The energy regime recently opened up by the advent 
of the LHC offers new possibilities. Clearly pA data 
at the LHC are essential to investigate various physics 
effects, in particular nuclear shadowing in a still unex- 
plored x-range. The final-state interaction of quarkonia 
with CNM is another attractive topic. Due to the ex- 
pected very short overlap time of the heavy-quark pair 
with nuclear matter, one could expect its influence on 
the still almost point-like cc to be rather hmited. How- 
ever, these considerations are still very quahtative and 
deeper theoretical studies are clearly necessary. 

The study of open heavy flavour production in pp 
and pA colUsions is tightly cormected to the process of 
quarkonia (hidden flavor) production. It offers a good 
basis to estabhsh fundamental parameters for quarko- 
nium calculations, together with an outstanding frame- 
work to test fundamental concepts, such as perturbative 
QCD, factorization and non-perturbative power correc- 
tions. Because heavy quarks are massive, the produc- 
tion cross section can be calculated in perturbative QCD 
down to Pt = at the parton level. 

However, differential distributions and event shapes 
exhibit large logarithmic contributions, typically cor- 
responding to soft or colhnear parton radiation, which 
need to be resummed to all orders. Furthermore, 
dead-cone effects suppress gluon radiation around the 
heavy-quark direction, which has relevant phenomeno- 
logical impUcations on both parton- and hadron-level 
spectra. As far as heavy-hadron production is con- 
cerned, one can model non perturbative effects by 
means of hadronization models such as the Lund-based 
string models, which are implemented within the frame- 
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work of Monte Carlo generators HERWIG, PYTHIA, 
MC@NLO, etc. Multi-purpose generators have been 
available for several years for lepton/hadron collisions 
in the vacuum and some have been lately modified to 
include the effects of dense matter. 

In order to promote the formalism used to describe 
open heavy flavors to pA and ultimately AA collisions, 
one has to introduce medium-modification effects. In 
particular, one of the most striking observations of 
heavy-ion collisions is the jet-quenching phenomenon, 
namely the suppression of hadron multiplicity at large 
transverse momentum with respect to pp processes. In 
the case of heavy hadrons, the measured suppression is 
expected to be lower than that of light-flavor mesons, 
due to the dead-cone effect. So far, suppression of 
heavy-flavored mesons in pA or AA collisions has only 
been measured through leptons (non-photonic electrons 
and single muons), i.e., the leptons from semi leptonic 
decays of heavy-flavored hadrons. Surprisingly, the en- 
ergy loss observed through heavy-flavor leptons is sub- 
stantial already at moderate pj of ~3 GeV/c. Several 
mechanisms to describe heavy-flavor energy loss in the 
dense medium have been proposed. Most prominent 
are coUisional and radiative energy loss, typically calcu- 
lated in the weak-coupling regime. However, the mag- 
nitude of the observed suppression at RHIC is hard to 
accommodate in these models. Recent attempts, assum- 
ing strong coupling between the heavy-flavor mesons 
and the dense medium and based on AdS/CFT corre- 
spondence lead to higher suppression factors. Further 
studies as well as experimental data are needed before 
any conclusions can be drawn. Future data from the 
LHC wiU possibly shed light on the issue of the energy 
loss. It will be very useful to have more exclusive probes 
such as azimuthal correlations or tagged c- or ^-flavored 
jets. 

Another crucial issue in heavy-flavor phenomenol- 
ogy concerns heavy-quark parton distributions in pro- 
tons and nuclei. In fact, processes with the production 
of heavy quarks are fundamental to constrain the PDFs. 
In particular, the production of a direct photon accom- 
panying a heavy quark at hadron colliders is an impor- 
tant process, also relevant to constrain the heavy-quark 
density. As it escapes confinement, a photon can be ex- 
ploited to tag the hard scattering; also, its transverse mo- 
mentum distribution gives meaningful information on 
the heavy-quark production process and PDF. 

The LO/NLO calculations for y + c/b production can 
be extended to pA collisions, e.g. proton-lead processes 
at the LHC. To a first approximation, the process can be 
calculated by convoluting the vacuum matrix elements 
with nuclear PDFs. While the gluon distribution in the 



proton is rather well constrained, it is poorly known 
in the nucleus, hence large discrepancy between gluon 
shadowing parametrizations. In fact, possible measure- 
ments at the LHC of y + c/b final states will help to 
investigate both gluon and charm/bottom distributions 
in dense matter. Since the Bjorken-x regions probed 
at RHIC and LHC are complementary, these measure- 
ments will help to discriminate among the nuclear par- 
ton distribution sets. 

Lattice QCD calculations predict that, at sufiiciently 
large energy densities, hadronic matter undergoes a 
phase transition to a plasma of deconfined quarks and 
gluons (QGP). Considerable efforts have been invested 
in the study of high-energy heavy-ion collisions to re- 
veal the existence of this phase transition and to study 
the properties of strongly interacting matter in the new 
phase, in view of e.g. improving our understanding of 
confinement, a crucial feature of QCD. The study of 
quarkonium production and suppression is among the 
most interesting investigations in this field. Calcula- 
tions indicate that the QCD with binding potential is 
screened in the QGP phase, the screening level increas- 
ing with energy density. Given the existence of several 
quarkonium states of different binding energies, it is ex- 
pected that they wiU be consecutively dissolved (into 
open charm or bottom mesons) above certain energy 
density thresholds. 

One of the objectives of high energy nuclear colli- 
sions is to produce and study the QGP under controlled 
conditions in the laboratory. How then can we probe 
the QGP - what phenomena provide us with informa- 
tion about its thermodynamic state? Here the ultimate 
aim must be to carry out ab initio calculations of the 
in-medium behaviour of the probe in finite tempera- 
ture QCD. Quarkonia may well provide the best tool 
presently known for this. However, quarkonium pro- 
duction wiU also be modified in nuclear collisions by 
other phenomena. Parton distribution modifications, 
parton energy loss and cold nuclear matter effects on 
quarkonium production must be accounted for before 
any QGP studies become meaningful. 

Any modifications observed when comparing 
J/i// production in AA collisions to that in pp interac- 
tions thus have two distinct origins. Of primary interest 
is obviously the effect of the medium produced in the 
collision - this is the candidate for the QGP we want 
to study. In addition, however, the presence of the 
"normal" initial and final state effects will presumably 
also affect the production process and the measured 
rates. This ambiguity in the origin of any observed 
J/>// suppression wiU thus have to be resolved. 

A second empirical feature to be noted is that the 
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observed Jjij/ production consists of directly produced 
\S states as well as of decay products from xAXP) and 
i/''(2S) production. The presence of a hot QGP affects 
the higher excited quarkonium states much sooner (at 
lower temperatures) than the ground states. This results 
in another ambiguity in observed Jjil/ production - are 
only the higher excited states affected, or do all states 
suffer? 

Ideally, one would measure Xc and ^' produc- 
tion separately first in pA (or dA) collisions, to deter- 
mine the effects of cold nuclear matter, and then, again 
separately, in AA collisions as a function of collision 
centrality at different center-of-mass energies. 

To deepen our understanding of heavy-quarkonium 
production mechanisms inpp and pA, we must now be- 
gin to study new observables such as associated produc- 
tion. As we have discussed above, theoretical uncertain- 
ties on the yields are too large to draw any strong con- 
clusions in favour or disfavour of one model or the other. 
In addition, polarization studies are in practice chal- 
lenging and not always possible in the complete phase- 
space region where one can measure the yield. Much is 
expected from the LHC on the latter two observables. 
Nonetheless, one should also use the LHC for exper- 
imental investigations of new measurements, not pos- 
sible previously. These observables can test the many 
production models available, in particular the CSM, the 
COM, the CEM and the kj factorization. The study of 
associated production, ijj + cc,l' + bb, ij/ + y, . . . , first in 
pp collisions, then in pA and AA, should be very fruit- 
ful. 
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